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Foreword 


In  an  attempt  to  stimulate  communication  among  the  many  groups 
vorV.ing  on  droplets,  sprays  and  atomization,  Project  SHU1D  undertook 
to  sponsor  a  conference  devoted  to  this  area  of  study.  The  confer¬ 
ence  took  place  at  Northwestern  University  on  September  21*  and  25, 

1955-  Most  of  the  work  involved  in  preparation  and  staging  was  ably 
carried  out  by  Professor  George  Lamb  and  his  associates  in  the  Chemical 
Engineering  Department  of  the  Northwestern  Technological  Institute. 

S'«UID  takes  this  opportunity  of  registering  its  appreciation  for  these 
efforts  ar.d  for  the  cooperation  and  hospitality  offered  by  the  University. 

As  initially  conceived,  the  conference  was  to  be  devoted  largely  to 
more  or  less  informal  review  and  discussion  of  various  programs  and  topics 
relating  the  formation  of  sprays  and  droplets,  their  properties,  and  their 
applications.  A  primary  objective  was  to  have  as  complete  representation 
as  possible  of  all  the  groups  active  in  spray  and  droplet  research  to¬ 
gether  with  representatives  of  sponsoring  agencies  in  order  that  future 
vork  in  the  field  might  be  planned  with  some  discretion,  and  minimum 
duplication. 

Originally,  publication  cf  the  proceedings  of  the  Meeting  was  not 
contemplated.  It  was  believed  that  more  complete  participation  might  be 
achieved  by  encouraging  informal  discussion  and  by  sparing  contributors 
the  chore  of  preparing  formal  manuscripts.  Since  the  meeting  there  have 
been  a  large  number  of  requests  for  some  sort  of  summary.  This  report 
is  in  response  to  these  requests.  It  Includes  abstracts  of  most  of  the 
talks  and  complete  papers  for  four  of  them.  In  some  cases  it  has  not  been 
possible  to  obtain  even  abstracts.  The  render  is  referred  to  the  authors 
in  these  instances,  it  is  hoped  that  this  cunjr.ury  of  the  conference  pro¬ 
ceedings,  even  though  incomplete,  will  serve  a  useful  purpose. 
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Abstract 

This  report  is  a  sunvnury  of  the  proceedings  of  a  conference  on 
Atomization,  Sprays,  and  Droplets  sponsored  by  S'XU1D  at  or thve stern 
Technological  Institute  in  September  1955-  Seventeen  papers  and  dis¬ 
cussions  vere  presented.  Abstracts  or  references  to  each  are  given. 
Four  of  the  papers  art-  giver,  in  entirety. 
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September  PA-o,  iSp.5 

Paper  Number  1 

REVIEW  OF  EXPERIMENTS  AND  EMPIRICAL  COR RELATIONS 
hr  DATING  TO  THE  PRODUCTION  OF  SPRAYS 

by  W.  R.  Marshall,  Jr. 

University  of  Wisconsin 

No  abstract  of  the  paper  as  presented  Is  available.  However,  the 
r  iterial  covered  will  appear  as  part  of  a  Cherni  .tl  Engineering  Progress 
v-nogr aph  by  Dr.  Marshall  on  Ato~4 zation  and  Spray  Drying.  This  mono¬ 
graph  will  be  published  by  the  American  Institute  of  Chemical  Engineers 

and  is  now  in  process. 
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Paper  Number  2 


REVIEW  OF  THEORETICAL  AND  MATHEMATICAL  ANALYSES  OF  THE 

PERFORMANCE  OF  ATOMIZING  NOZZLES 

> 


By  R.  Laster  and  M.  Doumas 
General  Foodif  Corp. 
Hoboken,  N.  J. 


*•’0  abstract  of  this  paper  is  available.  The  material  is  included 
in  a  paper  by  the  authors  appearing  in  Chemical  Engineering  Progress 
Volume  I49,  Number  10,  pages  518-526,  October  1955* 
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S  t^U  ID  CONFE.ilT  ICE  CN  ATOMIZATION’.  SPRATS,  AND  DROP  UTS 
Northwestern  Technological  Institute 
--September  24-Sp,  1953 

Paper  Number  j 

CHARACTERISTICS  OF  SPRAYS  AND  DROPLETS 

By  J.  Mason  Pilcher 
Battelle  Memorial  Institute 

The  cascade-impactcr  method  has  been  the  subject  of  experimental 
study  at  Battelle  during  the  past  2  V2  years.  A  review  of  this  work 
i 8  presented.  If.  this  device,  droplets  in  a  moving  air  stream  impact 
upon  a  slide  placed  in  their  path  provided  the  inertia  of  the  particles 
is  sufficient  to  overcome  the  drag  exerted  by  the  air  stream  that  must 
move  around  the  slide.  As  the  mist  is  drawn  through  the  various  stages 
of  the  impactor,  tne  speed  of  the  stream  and  therefore  the  efficiency  of 
impaction  increase  and  a  size  classification  of  droplets  is  accomplished. 
The  Battelle  instrument  consists  of  six  Jet  stages  and  a  seventh  filter 
stage  arranged  in  tandem  fashion.  After  this  instrument  is  once  cali¬ 
brated,  rapid  determinations  of  drop  size  distribution  may  be  made  by 
vising  gravimetric,  colorimetric,  fluorometric  or  radiochemical  methods 
of  analysis,  instead  of  the  tedious  procedure  of  microscopic  counting. 

The  impactor  is  most  suitable  in  the  size  range  of  about  1  to  50  microns. 

The  movement  of  an  Individual  droplet  issuing  from  a  spray  nozzle  is 
analyzed.  In  this  review,  it  Is  shown  how  the  motion  depends  upon  the 
initial  velocity  of  the  drop  and  the  drag  due  to  the  surrounding  air. 

Drag  coefficients  of  reasonable  accuracy  are  available  in  the  literature; 
however  the  analysis  of  the  motion  is  greatly  complicated  by  the  period  of 
acceleration  and  the  distortion  of  the  drop. 

Additional  experimental  data,  especially  high-speed  photography  of 
droplets  ejected  at  high  velocities  are  required  for  a  better  knowledge 
of  the  dynamic  behavior  of  droplets.  The  splnnlng-diok  atomizer  is 
suggested  as  an  excellent  tool  for  use  in  study! n.:  the  trajectory  and 
velocity  of  droplets.  (Kotfj  Ccmplote  inper  in  Appendix  1) 
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l’uper  Number  4 

CHOI’  SIZE  DISTKIBUTlOflS  AND  THEIR  AVERAGES 

By  Ralph  E.  Thomas 
Bat telle  Memorial  Institute 

In  the  current,  literature,  there  are  several  types  of  distributions 
which  arc  of  Importance  In  describing  l.ic  size  consist  of  sprays.  These 
dintr Ibuttonu  arise  naturally  In  the  course  of  experimental  Investigations 
employing  different  tec  rmlqucc .  The  specific  distributions  frequently  em¬ 
ployed  in  droplet  statistics  are:  the  normal  distribution,  the  Log-normal 
distribution,  the  distributions  of  Her. In -Hammier  and  Nukiyoma-Tanasawa. 
Having  chor.cn  a  particular  type  of  distribution  function,  the  constant 
parameters  may  be  computed  so  that  the  curve  gives  the  best  fit  to  the 
data.  Analytical  methods  should  supplement  the  graphical  methods  in  the 
final  determination  of  the  values  of  the  parameters.  The  customary  use 
of  relative  error  as  a  measure  of  the  difference  between  theoretical 
and  observed  results  Is  not  adequate  when  probability  considerations 
are  Involved,  Unless  the  systematic  errors  can  be  determined  and  their 
effects  eliminated,  the  fundamental  problem  of  droplet  statistics  cannot 
be  solved. 

(Note:  Complete  Taper  In  Appendix  2) 
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Pacer  Number  5 

REVIEW  OF  MASS  AND  ENERGY  TRANSFER  BETWEEN 
_ LIQUID  DROPLETS  AND  SURROUNDING  GASES.  _ _ 

By  Charles  C.  Graves  and  D.  VI.  Bahr 
National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio 

Close  agreement  tetveen  predicted  and  experimental  evaporation  rates 
of  single  drops  in  quiescent  surroundings  has  been  reported.  The  Langmuir 
equation  has  been  verified  for  the  evaporation  of  small  drops  in  room  air. 
The  Fuchs  equation  for  evaporation  of  small  drops  at  low  ambient  pressures 
has  been  confirmed  experimentally.  The  evaporation  of  single  drops  in 
high  temperature  surroundings  is  limited  to  the  special  case  of  single 
drops  burning  in  quiescent  atmospheres.  Good  agreement  between  theory 
and  experiment  was  obtained. 

Semi-empirical  correlations  for  steady  state  forced  convection  evap¬ 
oration  from  single  drops  of  pure  liquid  are  available  in  terms  of  cither 
the  heat  or  mass  transfer  equations.  Drop  surface  temperatures  may  be 
determined  for  simultaneous  solution  of  these  equations  or  from  a  correla¬ 
tion  in  terms  cf  the  liquid  boiling  point  for  the  higher  temperature 
evaporation  conditions.  Theoretical  heat  or  mass  transfer  coefficients 
for  spheres  are  in  approximate  agreement  with  experiment. 

Some  data  are  available  on  the  effect  of  both  dissolved  and  sus¬ 
pended  solids  on  the  evaporation  of  water  drops.  Ilo  suitable  data  are 
available  on  the  evaporation  of  wid'-  boiling  range  hydrocarbons. 

Calculations  have  indicated  that  the  unsteady  state  evaporation  may 
represent  an  appreciable  portion  of  the  total  evaporation  time.  However, 
little  experimental  work  has  been  done  on  this  [ihr.se  of  the  evaporation 

process . 

Theoretical  analyses  of  fuel  spray  evaporation  have  boon  made  with 
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tbe  assumption  of  aero  relative  drop-air  velocities.  Little  experi¬ 
mental  vorb  on  fuel  spray  evaporation  has  been  done.  To  date  no 
correlation  betveen  single  drop  and  total  spray  evaporation  has  been 
obtained.  ' 

(3ote:  Complete  Paper  in  Appendix  3) 
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Paper  Number  6 


FUEL  SPRAY  SPREADING  AND  EVAPORATION^  JN 
A  FLOWING  AIR  STREAM 

By  Donald  V? .  Bahr 

National  Advisory  Committee  for  Aeronautics 
Levis  Flight  Propulsion  Laboratory 
Cleveland,  Ohio 


The  evaporation  rate  and  degree  of  fuel  spreading  of  gasoline -type 
cprays  were  Investigated  over  the  range  of  inlet-air  conditions  common 
in  raa-Jet  and  tailpipe  engines.  This  study  uas  conducted  by  injecting 
isooctane  contrastreara  from  a  simple -orifice  nozzle  into  air  flovlng 
through  an  8-iach  diameter  duct.  Both  the  total  and  liquid  fuel  distri¬ 
butions  across  the  duct  were  determined  at  various  positions  downstream 
of  the  fuel  injector  by  sampling  measurements.  Some  similar  data  were 
also  obtained  from  teats  in  a  l6-inch  diameter  ram-jet  engine. 

The  influence  of  each  of  sis  experimental  parameters  vao  system¬ 
atically  determined,  us  follows:  The  measurements  were  tui.cn  at  axial 
distances  of  5,  10,  and  l8  inches  from  the  nozzle  over  the  ranges  of 
inlet-air  temperatures  und  velocities  of  540-960°  R  and  100-550  feet  per 
second,  respectively.  The  ambient  air  pressure  was  varied  from  18  to 
55  inches  of  mercury  absolute,  nozzle  diameters  of  0.024)  0.055,  and 
O.Cll  inches  and  three  fr  ’  injection  pressure  differentials,  25,  55, 
and  Q5  psig  were  utilized. 

Over  these  ranges  of  conditions,  the  evaporation  rute  and  the  degree 
of  fuel  spreading  were  correlated  in  terms  of  the  six  experimental  para¬ 
meters.  The  extremely  large  influence  of  ntr  temperature  on  the  spray 
vaporization  was  demonstrated  by  this  analysis.  The  expression  relating 
the  fuel  spreading  to  the  experimental  variables  indicated  that  the  spray 
residence  time  was  tha  most  important  consideration  in  determining  the 
degree  of  spreading. 
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The  results  of  this  investigation  are  limited  to  the  case  of  up¬ 
stream  injection  of  isooctane  from  simple  orifice  noiiies.  However, 
the  relationships  can  probably  be  extended  to  ar.y  gasoline -type  fuel 
without  serious  error.  The  fuel  spreading  results  are  also  limited 
since  these  data  were  obtained  with  turbulence  less  than  that  to  be 
expected  in  engines.  Therefore,  for  application  to  engine  design, 
these  results  represent  the  minimum  degree  of  fuel  spreading.  The 
major  contribution  of  this  study  was  to  provide  an  indication  of  the 
relative  importance  of  ee.ch  of  the  inlet-air  and  fuel-injection  para¬ 
meters  on  the  evaporation  end  mixing  of'  fuel  injected  into  high  vel¬ 
ocity  airstreams . 
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Paper  Number  7 

REVIEW  I :>oXnOMENTAT ION  AND  S.P1  rlODS  CK  r.'/'i  ERLMr.NTAi.  oiuDY 

By  J .  L.  York 

University  ;.f  Michigan 


No  abstract  of  this  paper  is  available.  Of  pn-ticolar  interest 
was  the  discussion  of  the  electronic  counting  probe  technique  for 
studying  sprays  which  has  been  developed  by  Dr.  York.  The  reader  is 
referred  to  the  author  for  further  information. 
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Paper  Number  8 

AN  APPARATUS  FOR  SUSPENDING  SMALL  DROPLETS  OF  LIQUID  IN  SPACE 

By  J.  A.  Bolt 
University  of  Michigan 

This  paper  describes  a  piece  of  apparatus  that  is  being  used  to 
suspend  single  droplets  of  liquid  in  space .  The  method,  unlike  all 
previous  methods,  does  not  depend  on  physical  contact  between  the  drops 
and  a  solid  filament  but  obtains  the  necessary  forces  for  holding  the 
drop  from.  a.  rising  column  of  air  and  a  standing  ultrasonic  field.  The 
air  column  supplies  the  necessary  drag  force  to  overcome  the  force  of 
gravity  while  the  ultrasonic  field  provides  both  lateral  and  vertical 
forces  to  hold  the  drop  in  a  fixed  position. 

To.  generate  the  air  column,  a  vertical  wind  tunnel  was  designed 
which  consists  of  a  divergent  section  to  decrease  the  velocity  of  the 
supply  air,  a  calming  section  to  reduce  turbulence,  and  n  convergent 
nozzle  to  supply  a  low  turbulence  column  of  air  to  the  ceramic  trans¬ 
ducer  described  below.  Auxiliary  apparatus  associated  with  the  wind 
tunnel  are  a  pressure  reducing  valve,  surge  tank,  air  filter,  air 
heater,  and  flow  orifice. 

The  ultrasonic  field  is  generated  by  a  cylindrical  tube  barium 
titauate  ceramic  transducer.  The  tube  is  four  inches  long  and  has  an 
inside  diameter  of  one  and  five-sights  inches.  Two  windows  are  cut  in 
the  wall  so  that  the  drop  can  be  illumineted  and  photographed.  Power 
for  the  transducer  is  supplied  by  a  forty-five  watt  amplifier  which  is 
fed  by  an  audio  oscillator  with  a  range  to  40,000  cps. 

In  operation  the  ceramic  transducer  io  placed  above  the  outlet 
of  the  wind  tunnel  nozzle,  and  a  lov  velocity  air  atrtum  is  started 
through  the  transducer.  The  oscillator  is  then  varied  until  the  de¬ 
sired  resonant  frequency  of  the  tube  is  reached,  «t  uoout  37,000  cps. 
This  is  detected  by  a  needle  vhich  is  suspended  in  the  tube.  Before 
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resonance  is  reached,  the  needle  shove  no  preferential  position  as  it  la 
moved  within  tne  tube,  but  at  resonance,  the  r.eer^e  is  held  i^alte  flretly 

along  or.c  of  the  nodes  In  the  ultrasonic  f  ic.i.  A  slight  deflection  uf 

the  string  used  to  augr.er.d  the  needle  will  act  ms-c  the  nccdlu. 

When  the  tube  is  at  proper  resonance,  droplets  of  liquid  ere  chases 
into  the  ceramic  tune  from  an  eye  dro j.per,  and  the  air  velocity  ic  ad- 

4  _  .... - S.  S  SU.  «-  ..4.  «..n4  #■«  aVaAvoe 

Uutil  a  tu  VV  w  ««4^Vi4Utu  S>|  OUV  u  ul  u  v-  <-  l>  iiyeu  eeSu  «  V .  VWaSWs  tta- 

tlon.  As  the  drop  evaporates,  the  air  velocity  is  r-anuaily  aecr  _»;£d 
to  keep  the  drop  In  position  opposite  the  obcervaticn  window. 

Photographic  records  of  evaporating  drops  are  made  with  an  Array 


frames  per  second.  To  obtain  the  necessary  magnification,  an  extension 
tube  for  the  lens  has  been  built,  giving  a  magnification  of  about  hX. 

The  extension  tube  also  contains  a  half  silvered  mirror  and  ground  glass 
so  that  the  drop  can  be  observed  as  it  is  being  photographed. 


HtASfRLMnKr  OF  Ml  bULE  TEXPESATURjE  IV  MISTS 

- -ry  *'«  t#  . -trcsinjs  - - 

i’rfiiio:  University 


An  analysis  of  the  evaporation  of  a  water  rpray  In  a  simplified 
flu-  which  b-sr  ir.  fact  a  prototype  of  »  ce«  design  >-*f  spray 

dryerj  is  presented.  In  addition,  a  new  instrument,  the  pneumatic 
thermometer  and  hygrometer  which  Is  useful  in  studies  ol*  sprays,  Is 
described.  Tbe  meaning  of  the  pressure  measured  by  a  Pltcfc  tube  in  a 
flowing  fcprdy-gas  mixture  is  analyzed  mathematically. 

The  water  spray  was  evaporated  In  a  prototype  of  the  high  velocity 
6pray  dryer.  This  dryer' consisted  of  a  tvo-fluld  atomizer,  vith  hot  air 
atomizing  the  liquid  feed,  vhich  discharged  into  a  second  concurrent 
stream  of  hot  air.  It  represented  an  attempt  to  utilize  the  rapid  heat 
and  mass  transfer  during  atoaizatlon  in  spray  drying.  The  momentum, 
mass  and  heat  transfer  in  this  system  were  followed  in  detail  by  measur¬ 
ing  the  velocity,  gas  temperature,  vater  temperature,  liquid  water  con* 
csntration,  ana  humidity  at.  a  large  number  of  points  in  each  of  six 
successive  cross-sections.  The  results  of  these  measurement  for  a  single 
case  of  water  evaporating  Into  air  are  presented  and  analyzed.  Water  was 
atomized  to  droplets  cf  20  microns  moan  diameter  by  a  stream  of  air  at 
l*00°f ' .  moving  at  8;0  ft  ./sec.  This  stream  was  subsequently  mixed  with  a 
stream  at  280°F .  moving  at  50  ft. /sec.  The  atomizing  nozzle  was  0-500 
inches  in  diameter  and  tbe  coaxial  drying  chamber  was  6  inches  in  Jia-  ■ 
meter.  Evaporation  of  the  spray  was  complete  within  50  inches  down¬ 
stream  from  tbe  nozzle.  Tbe  data  were  analyzed  in  terms  of  Relchardt's 
theory  of  free  turbulence.  The  spreading  coefficients  obtained  were  in 
the  same  range  as  those  obtained  by  other  investigators,  although  not  in, 
close  agreement.  The  spreading  coefficient  for  water  droplets  was  measured 

-1/t- 


and  found  to  be  considerably  dess  than  that  for  gas  under  the  same 
conditions . 

The  pneumatic  thermometer  and  hygrometer  was  developed  to  permit 
the  measurement  of  the  gas  temperature  end  humidity  in  a  gas-spray  mixture 
This  instrument  relies  on  the  changes  in  pressure  and  temperature  as  a 
sample  of  tht  mixture  is  passed  through  as  orifice,  then  a  heater,  a 
second  orifice,  a  drying  bed,  and  finally  a  third  orifice.  Details  of 
the  development  and  testing  of  .nis  new  instrument  are  presented. 

A  mathematical  analysis  of  the  performance  of  a  Pitot  tube  in  a 
flowing  spray-gas  mixture  is  presented.  This  analysis  gives  a  method  of 
determining  the  magnitude  or  tne  error  caused  Dy  tne  droplets  entering 
the  Pitot  tube.  It  ie  shown  that  this  error  is  dependent  upon  drop  sice 
and  Pitot-tube  diameter,  and  that  a  series  of  Pitot  tubes  may  therefore  be 
used  to  determine  drop  size  distribution  in  a  gas-boms  spray.  Su^h  a 
device  awaits  practical  development. 

Microfilm  copies  of  the  thesis  are  available  through  University 
Microfilm,  Inc.,  Ann  Arbor,  Michigan. 


Note:  The  material  covered  by  Dr.  Comings  was  taken  from  a  thesis 
entitled  "Sprays  in  Hot  Turbulent  Streams"  by  Clarke  Lincoln 
Coldren,  Department-  of  Chemistry  and  Chemical  Engineering, 
University  of  Illinois.  The  thesis  lias  been  deposited  with 
the  University  of  Illinois  Library.  Microfilm  copies  are 
available  through  university  Micr->filrc  Inc..,  Ann  Arbor,  Michigan. 
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Paper  Kicrbcr  10  i 

CAS  ABSORPTION  AND  AERCSoI  '  COLLECTION  IN  A  VENTURI  ATOMIZER  | 

By  H.  F.  Johnstone,  R.  B.  Feild1,  and  M.  C.  Tassler^ 

University  of  Illinois 
Urbar.a,  Illinois 

I 

The  success  of  the  Venturi  scrubber  as  an  efficient  device  for 
contacting  large  volumes  of  cases  for  absorption  and  dust  collection  ! 

*  •  >  I 

has  made  important  this  study  of  the., rates  and  mechanism  of  transfer 

in  the  spray  zone  .  The  scrubber  consists  of  a  conventional  Venturi 

section  into  which  the  scrubbing  liquid  is  injected  at,  the  throat, 

followed  by  a  cyclone  for  removing  the  liquid  droplets  containing  the 

material  collected.  The  high  velocity  of  the  gas  stream  relative  to  the 

liquid  causes  atomization  into  a  spray  of  many  small  droplets  which  present 

a  large  surface' for  mesa  transfer.  | 

i 

Gas  film  coefficients  were  measured  for  the  absorption  of  sulfur  1 

dioxide  by  t n  alkaline  solution  by  collecting  the  spray  droplets  by  an 

impact  nozzle  surrounded  when  necessary  by  a  stream  of  inert  gas.  Ex-  ! 

posure  tiir.es  for  the  liquid  ver.e  in  she  range  from  0.1  millisecond  to  i 

1  secund.  Similar  methods  were  used  for  measuring  the  liquid  film  co¬ 
efficients  for  the  desorption  of  carbon  dioxide  from  vatec,  and  the  rate 
of  collection  of  aerosol  particles  by  water  spray.  The  aerosols  were  pro¬ 
duced  by  the  condensation  of  dibuty Iplithalate  vapor,  or  by  the  chemical 
reaction  of  ammonia  gas  with  dilute  SQ^  or  HC1  in  air.  The  particle  size 
of  the  aerosols  were  nearly  uniform  and  varied  between  0.27  and. 1.22  micron 
mean  diameter  for  the  three  materials  used. 


^Present  address:  Eastern  Laboratory,  du  Pont  Company,  Glbbstoun,  New  Jersey. 
“Present  address:  Film  Division,  du  Pont  Cor-pany,  Buffalo,  New  York. 
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The  Venturi  section*  used  were  made  cf  Lucite  with  throats  1-1/8  to 
1-1/2  inchej  in  diameter  and  1-1/2  to  j  inches  long  and  vjth  U°  to  70, 
divergence  angles . 

Che  observed  individual  absorption  film  coefficients,  kLa  and  k^-a, 
are  very  high  near  the  liquid  scarce.  The  value  oi  the  former  for  COq 
falls  off  abruptly  from  a  maximum  about  one  inch  downstream  which  is  in 
the  range  of  5000  to  17;.OQ0  lb.  molej7'(hr.)(cu.xt.)(lb.  mcie/cu.rt .)  to 
almost  zero  at  four  inches.  The  average  gas_ film  coefficient  for  SO2  is 
in  the  range  o£_3D0  to  lbOO  lb.  molcs/(hr.)(cu.ft.)(atm.)  at  3-5  inches 
from  the  liquid  source  and  decreases  more  slowly  vith  distance  as  the 
interfacial  Velocity  decreases.  Both  coefficients  are  functions  of  gas 
velocity  and  liquid-gas  ratio.  Estimates  of  the  surface  area  of  the 
liquid  droplets  based  on  atomization  equations  permit  calculating  the 
film  coefficients  kr  and  kg,  and  indicate  the  extremely  rapid  transfer 
rates  in  the  atomization  zone.  Possible  applications  of  Venturi  atomizers 
are  suggested  for  difficult  absorption  problems. 

Aerosol  collection  by  spray  droplets  In  the  atomization  zone  takes 
place  primarily  by  the  inertial  impaction  mechanism.  The  efficiency  of 
collection  depends  on  the  Inertia  of  the  aerosol  particles  and  is  pract¬ 
ically  independent  of  small  natural  electrostatic  charges  on  the  droplets 
and  the  particles,  and  of  their  diffusivity  and  wettability.  The  rate  of 
collection  falls  off  rapidly  with  distance  from  the  spray  nozzle  as  the 
relative  velocity  of  the  gas  and  liquid  decreases.  A  theoretical  equation 
based  on  inertial  for.ces  gives  excellent  correlation  of  the  collection 
efficiencies  for  all  of  the  aerosols  studied. 
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Paper  Number  11 


FUEL  INJECTION  PROBLEMS  IN  TURBOJETS 


By  H .  A  r  Fremont 
General  Electric  Company 
Aircraft  Gas  Turbine  Division 
Cincinnati  15,  Ohio 


No  abstract  of  this  paper  is  available.  Mr.  Fremont  gave  a  talk, 
illustrated  by  slidcc,  discussing  in  general  the  technology  of  turbo¬ 
jet  fuel  combustor  applications,  a  feu  of  the  present  problene,  and 
coma  areas  in  which  future  research  is  needed. 
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Paper  Humber  12 

INJECTION  IN  AFTERBURNING  TURBOJETS 

By  D.  G.  Samar&a  7  . - 

Wright  Patterson  A.ir  Force  Base 
Dayton,  Ohio 


1.  The  injection  of  fuel  in  afterburning  turbojets  has  a  profound 
effect  upon  the  performance  and  the  operational  stability  of  the  after¬ 
burner. 

2.  Generally  speaking,  at  least  two  operational  cases  of  after¬ 
burners  are  apparent: 

a.  The  afterburner  operates  only  a  short  time  of  the  total 
flight  time. 

h.  The  afterburner  operates  during  most  of  the  whole  flight 
tine.  ' 

3.  The  design  requirements  of  the  above  distinct  cases  are  entirely 
different. 

k.  The  types  and  geometric  locution  of  the  injectors  in  the  after¬ 
burner  have  a  great  Influence  upon: 

,  a.  Combustion  efficiency 

b.  Pressure  losses 

c.  Temperature  distribution 

d.  Instability 

5«  The  two  main  -lodes  of  upstream  and  downstream  injection  may  re¬ 
sult  in  the  following  five  types  of  afterburning: 

a.  Preturbine  injection 

b.  Afterburning  upstream  injection  with  flamehoiders 

c.  Afterburning  upstream  injection  without  flamehclders  (turbine 
stabilization) 


d.  Afterburning  downstream  injection  with  flams  holders 

e.  Afterburning  downstream  injection  with  tail-cone-end 
injector 

6.  The  advantages  and  disadvantages  of  the  above  types  are  disclosed 
A  hypothesis  of  the  screech  instability  is  given,  and  the  influence  of  the 
injection  process  upon  this  instability  is  discussed.  The  influence  of 
injection  during  acceleration  and  deceleration  <s  also  dealt  with. 


Note:  Complete  Paper  in  Appendix  4 
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Paper  Number.  13 

FUEL  INJECTION  PROBLEMS  IN  RAMJETS 

By  J,  P.  Longwell  „ 

Standard  Oil  Development  B 

Ramjet  fuel  is  typically  injected  into  the  high  velocity  air  at  the 
diffuser  exit  through  an  atemizing-r.ozzle  of  coderate  pressure  drop.  The 
air  stream  velocity  Is  generally  several  hundred  ft. /sec.  so  tudt  tuc  ic" 
lative  velocity  between  fuel  and  air  and  therefore  the  atomization  i3 
largely  controlled  by  the  air  velocity.  Air  atomization  data  using  pres¬ 
sure  atomizing  nozzles  for  a  variety  of  air  densities  are  therefore  of 
particular  interest.  The  drops  formed  are  in  the  p-LO  micron  range.  Such 
drops  are  capable  of  following  some  of  the  turbulent  velocity  fluctuations 
and  accelerate  very  quickly  to  stream  velocity.  The  rate  of  evaporation 
and  turbulent  diffusion  of  such  small  drops  requires  further  study  and  may 
not  he  easily  predicted  from  data  obtained  cn  large  drops. 
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'  APPLICATION  OF  SPRAY  TECHNOLOGY  TO  LIQUID  PROPELLANT  ROCKETS 

By  C.  C.  Mlccse 
Aerojet-General  Corporation 

Azusa,  California  - — 


Combustion  in  the  thrust  chamber  of  a  liquid  propellant  rccket  is  a 
very  complex  phenomenon,  so  that  it  is  almost  impossible  to  determine 
which  oi  the  many  processes  which  occur  is  the  determining  factor  for 
performance,  heat  transfer,  or  stability.  However,  on  the  basis  of  the.  1 
time  elapsed  during  each  of  the  processes,  it  is  logical  to  ess’imo  that 
the  less  rapid  physical  processes  of  jet  disintegration,  evaporation  and 
diffusion  should  have  c  greater  effect  than  the.  rapid  chemical  conversion 
of  the  unburned  propellants  into  the  products  of .. combustion.  It  la  inter¬ 
esting  to  note  that  the  theories  of  combustion  stability  -  us  presented  by 
Gunder  and  Friant,  Summer field,  and  Crocco,  for  instance  -  are  all  based, 
on  the  concept  of  the  time  lug  which  exists, between  tbs  injection  of  the 
propellants  and  their  subsequent  conversion,  thus  giving,  substantial  con¬ 
firmation  to  the  above  assumption  regarding  the  relative  Importance. of  the 
physical  processes. 

:  A-  •  -  •'  • 

Recognizing  the,  importance  of  the  physical  processes,  the  Aerojet-  , 
General  Corporation  has  undertaken  3eve.rul  experimental  and  theoretical 
investigations  of  Jet  disintegration,  droplet  evaporation,  and  gaseous 
diffusion. 

An  experimental  and  theoretical  study  of  jet  disintegration  in  ambient 
air  resulted  in  soire  interesting  correlations,  on  the  basis  of  previously 
published  theories,  of  data  obtained  from  two  types  of  injectors  (similar 
to  conventional  rocket  injectors),  and  using  two  different  liquids  -  water 
and  liquid  nitrogen.  The  latter  liquid  was  chosen  in  order  to  6tudy  vhat 
effects,  if  any,  were  caused  by  a  high  evaporation  rate,  and  it  was  found 


that,  the  apparent  dif ferences  between  the  water  and  liquid  nitrogen  data 
CQU-d  o e  ascribed  solely  to  the  dl fference  between  primary  and  disinte¬ 
gration  and  secondary  atom; tat  ion.  Th-  data  included  wave-length  of  dis¬ 
turbance,  maximum  drop-size,  and  break-up  length  of  the  Jet,  Additional 
experiments  are  currently  under  way  to  detersilnf;  the  mutual  affects  of 
several’  a oil near  j*ta, .and  the  effect  of  variations  in  ambient  pressure. 


A  second  phase  of  the  study  has  been  a  theoretical  investigati  on  »f 
the  ballistics  of  burning  droplefg.  By  assuming  that  the  drag  coefficient 
varies  inversely  as  the  Reynolds  number,  that  the  velocity  of  the  surround*, 
ing  gases  varies  directly  with  distance  downstream  of  the  injector,  and  the 

Al1  *  Is  r»  4  »iv  d«aa  1  ng  ^  I  ^  v-  1  ,  *  *,  (  t  „  f  4  n>n  W(1H+  I  a  *!♦%»)■ 

lt«4S  4  WW'<  Us  i.u  W1  slit>  u  ui  ililif,  SSI  I/W  4W  w  •  V>  S  A ».  o  .aauvt**  a*w<4  «*wsj  Iims  wvu  w 

of  notion  was  used  to  determine  the  velocity  profile,  lifetime,  and  proper 
initial  size  of  the  droplets  for  various  initial  velocities  and  ratios  of 
kinematic  viocosity  to  evaporation  rate.  Application  of  the  some  general 
method  tc  the  similar  pi  obloics  of  evaporation  or  burning  in  a  gas  stream 
whose  velocity  varies  linearly  with  distance  lb  outlined,  / 


Two  aspects  of  the  effect  of  diffusion  on  the  stability  of  cL  rocket 
uysteui  are  Indicated,  and  it  is  concluded  that  the  needs  and  opportunities 
for  further  research  ore  essentially  unlimited. 
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VAPORIZATION  OF  FUEL  DROPLETS  IN  A  DIESEL  ENGINE 

2y  ?.  S.  Myers,  0.  A.  Uyehara  and  M.  M.  El  Vakil 
University  of  Wisconsin 
Madison,  Wisconsin 

*  "  Theoretic  al  relationships  were  developed  from  basic  heat -'transfer 

and  diffusion  theories  ov  which  the  history  of  a  fuel  droplet  injected 
inheated  air  can  be  predicted.  The  droplets  are  assumed  tu  be  injected 
with  a  certain  initial  velocity  relative  to  the  air  and  are  slowed  down 
by  aerodynamic  'drag  forcen,  the  .drag  coefficients  being  taken  as  those  for 
solid  spheres.  The  vaporization  time  is  roughly  divided  into  two  portione, 
the  unsteady  state,  or  the  time  spent  by  the-- droplets  in  heating-up  to  their 
'vet-bulb  temperatures ,  corresponding  to  the  air  ambient  conditions,  and  the 
steady  utatc  where  the  droplets  remain  at  the  vet-bulb  temperatures  until 
comp.late  vaporization  takes  place.  The  unsteady  state  portion  has  been  found 
to  be  of  importance  in  roost  cases. 

The  equations  assume  (a)  no  interaction  between  droplets  (b)  the  air 
ambient  conditions  remain  unchanged  throughout  the  vaporization  procesr.  and 
(c)  spherical  symmetry  uf  the  droplets  and  the  air -vapor  film  surrounding 
them  exists  at  ala  times.  The  equations  are  integrated  by  a  stepwise  pro- 

i  . 

endure  in  which  their-  lifetime  is  divided  into  small  Inc:  e:r.:nts  during 
which  the  droplet  radius,  temperature  and  velocity  relative  to  the  air  are 
assumed  constant.  It  was  found  that  the  equations  yield  rensonuciy  accurate 
results  When  compared  with  experimental  values  obtained  for  fairly  large 
drops  vaporizing  in  heated  atmospheric  air. 

While  the  droplets  are  vaporizing  i.n  the  unsteady  state,  they  reach 
such  a  temperature  where  their  vapor  pressure  is  high  enough  oo  that  a 
stoichiometric  mixture  of  air  and  vapor  is  first  forme a  in  the  film  at  the 
droplet  surface  and  having  a  temperature  equal  to  that  of  the  liquid. 


As  the  droplet  it  heated  more,  its  vapor  pressure  increases  and  the 
stoichiometric  mixture  moves  away  from  the  droplet  in  .0  a  higher  temper¬ 
ature  region.  The  temperature  of  the  stoichiometric  mixture  reaches  a 
steady  state  at  essentially  the  same  tine  the  droplet  reaches  i+s  uet- 
bul'o  temperature.  Thus,  the  stoichiometric  mixture  must  reach  its  ig¬ 
nition  temperature  chile  the  droplet  is  in  the  unsteady  state  if  the 

—  drop-lxt-  is  lo  ignite  at  all .  The  time  spent  by  the-draplet-  until  the—' - 

stoichiometric  mixture  reaches  its  ignition  temperature  is  defined  as 
the  vaporization  portion  of  the  total  physical  ignition  delay  period. 

This  portion,  d,  of  the  physical  ignition  delay,  calculated  by  the 
above  method  was  found  to  approximately  follow  an  equation  of  the  form 

d  “  e  T 

where  T  is  the  absolute  air  temperature  and  c2  is  the  slope  of  this  line 
vhen  plotted  on  semi-log  graph  paper,  c-,  was  found  to  be  in  the.  range 
of  the  slopes  of  the  physical  ignition  delays  obtained  graphically  from 
experimental  combustion  bemb  data.  It  is  hoped  that  continuation  of  this 
theoretical  study  together  with  experimental  bomb  data  on  pure  hydrocarbon 
fuel6  will  help  In  formulating  a  rational  theory  for  physical  ignition 
delays . 

In  developing  the  above  equations  the  thermal  conductivity  of  toe 
liquid  droplets  was  assumed  to  be  infinite.  This  was  verified  experi¬ 
mentally  by  the  absence  of  a  temperature  gradient  inside  the  drops.  Ex¬ 
planation  was  found  inthe  internal  circulation  within  the  drops  which 
were  observed  during  the  experiments.  This  circulation  induced  enough 
fluid  mixing  to  effectively  make  the  thermal  conductivity  infinite. 
Circulation  was  found  to  be  in  a  doughnut  shaped  pattern  with  the  fluid 
near  the  surface  of  the  drop  moving  in  the  sate  direction  as  that  of  the 
air  stream  and  in  the  opposite  direction  at  the  center.  A  film  showing 
the  circulation  was  prepared  under  ar.  XACA  conLract  and  is  or.  file  at  the 
KACA  Washington  Headquarters. 
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Paper  Number  l6 


FUEL  INJECTION  PROBLEMS  IN  OPTO  CYCLE  ENGINES 

By  VI.  I.  Meyers 
Pennsylvania  State  College 
State  College,  Pennsylvania 

No  abstract  of  this  paper  is  available.  Professor  Meyers  has 
kindry  agreed,  to  make  available  to  all  those  attending  the  conference 
who  request  it  a  copy  of  the  Bibliography  of  Sprays  which  was  sponsored 
by  the  Texas  Company. 
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SQUID  CONFERENCE  OS  ATOMIZATION,  SPRAYS,  AND  DROPLETS 
Northwestern  Technological  Institute 
September  2b-25,  1955 

Paper  Humber  17 

CHEMICAL  ENGINEERING  APPLICATIONS 

By  H.  P,  Johnstone,  W.  Rl  Marshall, 

J.  L.  York,  S.  W.  Ceilings  and  G.  G.  Lamb 

This  ltea  on  the  program  was  not  a  formal  paper  but  comprised  a 
round  table  discussion  centering  largely  around  the  question  of  in¬ 
ternal  circulation  within  droplets.  In  addition,  Professor  Lamb  re¬ 
ported  on  his  work  at  Northwestern  dealing  with  the  evaporation  of 
two  component  fluids  atomised  in  pressure  swirl  nozzles  and  injected 
into  stagnant  heated  air. 
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by 


J.  Mason  Pilclier 


Bsttelie  Memorial  Institute 
Guptcuibcr  19a5 


CHARACTERISTICS  OF  SPRATS 
AKD  DROPLETS 


tiy 

J,  Kaaon  Pilcher 

_ The  primary  purpose  of  atomizing  a  liquid  is  to  increase  it3 

surface  area.  This  Is  especially  true  when  the  licrold  to  be  atomized  is  a 
fuel  that  must  be  vaporized  for  intimate  mixing  with  air  before  combustion 
can  take  place*  The  most  important  physical  characteristic  of  a  spray, 
therefore,  i3  drop-size  distribution,  because  this  characteristic  determines 
surface  area,  rate  of  evaporation,  and  rate  of  combustion. 

The  extent  to  which  atomization  increases  surface  area  is  ••'ell 
illustrated  by  J,  R.  Joyce^  by  means  of  a  simple  example.  A  cubic  centimeter 
of  liquid,  in  the  shape  of  a  sphere)  would  have  a  diameter  of  1.2)*  cm,  and 
its  surface  area  would  be  l*.  83  sq  cm.  If  this  sphere  is  then  atcmiced  into 
one  million  equal-sized  small  spheres,  each  having  a  volume  of  10  cc,  each 
droplet  would  have  a  diameter  of  0.012!*  cm,  and  the  total  surface  are*  of 
the  one  million  drcplets  would  bo  1*83  sq  cm,  or  one  hundred  times  that  of 
the  original  single  sphere.  The  surface  area  per  unit  volume  is,  therefore, 
inversely  proportional  tc  the  diameter  of  the  droplets,  for  by  reducing  the 
drop  size  by  a  factor  of  ICC,  the  total  surface  area  is  increased  by  a 
corresponding  factor  of  100. 

Figure  1  shows  these  relations  schematically.  The  usual  methods 
of  atomisation  result  in  the  licuid  being  broken  up  even  mere  finely  than 
is  indicated  by  this  example,  and,  of  course,  a  wide  range  cf  drop  sizes 
results,  A  typical  fuel  spray  would  consist  of  droplets  ranging  in  sizes 

from  1  micron  (0.CC01  cm)  to  SCO  microns  (0»C5  cm),  and  each  cubic  centimeter 
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BASIS:  one  cubic  centimeter  of  LIQUID 


BEFORE  ATOMIZATION 

AFTER  ATOMIZATION 

o 

NUMBER  OF  DROPS  *  1 
DIAMETER  =  1.24  CM 

AREA  =  4.83  SQ  CM 

VOL.  PER  DROP  =  1  CU  CM 

•  #  •  •  f<  , 

•  /•* 

,  •  ,  ,  •  i  •  v 

•  *  )  .  •  '  4  .•  «  *  '•  • 

:  ••  •> 

NUMBER  =  1,000,000 

DIAMETER  *  0.0124  CM 

AREA  =  483  SQ  CM 

VOL.  PER  DROP  =  IxlO-6  CU  CM 

FIGURE  I.  EFFECT  OF  ATOMIZATION  ON  SIZE  AND  SURFACE 
AREA  OF  DROFLETS 


of  fuel  would  toe  divided  into  not  Just  one  nil  lion  tut  tens  of  trillions  of 
droplets  having  a  surface  area  of  over  a  thousand  aq  cm  per  ce. 

The  instruments  and  experimental  methods  for  measuring  this  im¬ 
portant  characteristic  of  drop  size  will  be  covered  in  a  review  paper 
tomorrow  morning  t>y  Fruf  esso?  J.  1,  Yorkj  hf/wSV  5?  j  I  shall  include  a  brief 
discussion  of  only  one  method,  the  ciscade-itupactor  net  hod,  which  ha3  been 
the  subject  of  experimental  study  at  Battalia  during  the  past  2-1/2  years. 
Immediately  following  this  paper  Ralph  E.  Thomas  will  discuss  the  rather 
complex  and  controversial  subject  of  mathematical  representations  of  drop- 
size  distributions. 

The  other  major  characteristic  of  sprays,  in  addition  to  drop 
size,  is  the  dynamic  behavior  of  the  droplets  after  they  leave  the  atomizing 
device.  The  characteristic  of  dynanie  behavior  involves  such  factors  as 
the  initial  and  terminal  velocity  of  individual  particles,  the  time  and  dis¬ 
tance  of  travel  before  terminal  velocity  is  attained,  and  the  effects  of 
particles  upon  each  other  in  a  concentrated  spray. 

First  the  experimental  work  with  the  cascade  impactor  will  be  dis¬ 
cussed  briefly,  after  which  the  dynanie  behavior  of  droplets  will  be  reviewed 


MEASUREMENT  OF  DRPP-SI2E  DISTRIBUTION 
V.TTH  TPE  CASCADE  IMPACTOR  ' 

As  a  part  of  the  program  of  research  on  the  burning  characteristics 
of  fuel  mists,  extensive  work  has  been  done  on  the  development  of  the  cascade 
impactor  as  an  instrument  for  determining  the  drop-size  distribution  of 
atomized  liquid  fuel.  The  casoade-impaotor  method  was  chosen  from  among  the 
various  methods  that  you  will  hear  described  tomorrow  morning  because,  after 
the  instrument  is  once  calibrated,  rapid  determinations  of  drop-size 


distribution  nay  be  made  by  using  gravimetric,  colorimetric,  fluorometric, 
or  radiochssical  methods  of  analysis,  instead  of  following  tie  tedious  ard 
tir.e  consuming  procedure  of  microscopic  counting. 

Principle  Cf  The  Cascade  Tnpactor 

"'ha  cascade  inpactor  operates  on  the  principle  that  droplets  in  a 
moving  air  stream  rill  impact  upon  a  slide  placed  in  their  path  provided 
the  inertia  of  the  particles  is  sufficient  to  overcome  the  drag  exerted  by 
the  air  stream  that  must  move  around  the  slide* 

Figure  2  shows,  achecatically,  what  happens  when  a  fine  spray 
approaches  a  glass  slide  in  one  of  the  early  stages  of  a  cascade  lrpactor. 
The  large  droplets  impact  on  the  slide,  whereas  the  small  droplets  follow 
the  air  stream  into  the  next  Jet  which  is  smaller  in  diameter  than  the  pre¬ 
vious  one  and  which  thereby  imparts  a  higher  velocity  to  the  stream  of  mist 
At  the  higher  velocity,  smaller  sices  of  droplets  will  impact  on  the  slide. 
The  cascade  impactcr,  therefore,  acccm  oils  hen  a  size  classification  of  drop 
lets  by  increasing  the  speed  and  therefore  the  efficiency  of  impaction  as 
the  mist  is  drawn  through  the  various  stages  of  the  impactor,  K,  15.  Kay^, 
in  his  pioneer  work  on  the  cascade  impactor,  used  an  instrument  consisting 
of  four  stages  placed  at  right  anrles  to  each  other. 

The-  I-attelle  ?'o.  b  Cascade  Impactor 

The  instrument  now  being  used  at  “attolle,  by  J.taes  1.  Harp,  con¬ 
sists  of  six  jet  stages  and  a  seventh  filter  stage,  or  raved  ir.  tun  lorn  i'ash 
A  v— r*dcii  nry?. r.™?r.gnt  y*as  Tirst  \is?0.  by  t ,  t. 'r.n  ¥ 
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FUEL  MIST 


LARGE 


FIGURE  2. 


SCHEMATIC  DIAGRAM  SHOWING  PRINCIPLE 
OF  THE  CASCADE  IMPACTOR 
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Figure  3  shons  an  assembly  drawing  of  this  instrument  which  is 
designed  to  cover  the  size  range  of  3  to  1*0  microns  by  impacting  the  droplets 
on  the  six  Jet  stages.  A  seventh  filter  stage  removes  most  of  the  minus 
j -micron  material,  a  vacuum  applied  at  the  exit  of  the  irapactor  draws  the 
sample  through  the  instrument  at  a  fixed  rate  determined  by  a  critical  ori¬ 
fice  meter. 

Sampling  of  the  spray  is  a  difficult  problem  when  us:.*g  the  cascade 
impactor  and,  for  that  matter,  when  using  any  apparatus  for  determining 
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representative  as  possible,  an  "isokinetic"  sampling  procedure  is  followed. 
Isokivistic  sampling  is  achieved  when  the  velocity  of  the  mist  entering  the 
sampling  port  of  the  impactor  is  maintained  essentially  eaual  to  the  velocity 
of  the  stream  being  sampled.  The  path  of  the  mist  remains  practically 
unchanged  ss  it  enters  the  impactor,  and  die  crimination  against  eithsr  the 
larger  or  smaller  droplets  Is  reduced  to  a  minimum.  A  stationary  mist  ob¬ 
viously  cannot  be  sampled  lsokinetically. 


Impactor 


The  objective  in  calibrating  the  cascade  impactor  is  to  establish 
the  upper  and  lower  size  limits  or  3ize  range  of  the  droplets  that  collect 
on  each  stage.  This  is  accomplished  by  sampling  a  mist  over  a  short  enough 
period  of  time  to  prevent  coalescence  of  droplets,  and  measuring  the  diameter 
of  all  droplet  on  a  slide  or  portion  of  a  slide,  i’hotomicrographs  of  the 
impaction  slides  may  be  projected  by  means  of  an  enlarge  r-vicer  to  simplify 
greatly  and  improve  the  accuracy  of  calibration,  which  involves  the  measure¬ 
ment  of  about  a  thousand  droplets  Fer  Stage.  A  histogram  for  each  slide  is 


tiiea  prepared  to  establish  the  "cut-off1'  point  which  is  the  lerer  3ice  Halt 
for  one  stage  and  the  upper  size  Unit  for  the  succeeding  stages 

Cue  of  the  caj-'r  problems  in  the  calibration  of  the  iapacior  is 
that  the  cut-off  point  is  r.ct  perfectly  sharp.  end  a  certain  a:acur.t  of  over¬ 
lap  of  drop  sizes  occurs  between  adjacent  stages.  This  problem  is  resolved 
by  defining  trv*  e*tt-cfi‘  point  as  the  sice  of  drop  such  that  an  equal  nurd. or 
of  them  Mill  deposit  on  two  successive  stages. 

Other  factors  involved  in  calibrating  the  cascade  inpactor  arei 
(1)  the  use  of  a  nonvolatile  liquid,  3uch  as  dibutyl  phthalate,  and  (2)  cor¬ 
rection  for  the  flattening  of  spherical  droplets  after  impaction  on  the 
glass  slide.  Allowance  for  the  flattening  of  the  droplets  la  wo  do  by  de- 
termining  the  "flattening  coefficient",  which  is  the  factor  by  which  the 
flactened  diameter  mist  be  multiplied  to  obtain  the  diameter  of  the  spherical 
drop  before  Impact.  The  value  of  the  flattening  coefficient  for  droplets 
of  dibutyl  phcbalate  on  clean  glass  slides,  treated  with  a  retting  agent, 
is  0.?0  1  0.01  based  on  determinations  by  three  independent  ..methods. 

Once  the  cascade  iupactcr  is  calibrated,  and  the  upper  and  lover 
size  limits  for  each  slide  are  established,  rapid  determinations  of  drop- 
size  distribution  say  be  made  ty  employing  relatively  rapid  amlytical  methods 
to  determine  the  weight  of  material  impacted  on  each  stage,  of  a  tracer 

method  *  nonvolatile  dye  os*  fluorescent  22-lt  is  re c emended 

the  eternized  liquid  is  appreciably  volatile  as  is  the  case  for  all  liquid 
fuels. 


A  Limitation  Of  The  Cascade  Impact or 


is  that  it  docs  not  function  satisfactorily  when  sampling  coarse  sprays. 
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She  imp  Jeter  Is  -cat  suitable  in  the  site  range  of  about  1  to  SO  microns. 

T ha  high  inertia  of  large  particles  of  ICO  fiicrcns' diameter  ray  cause  them 
to  deviate  markedly  friz;  the  ?ir  flor  line;  before  the  first  silts  is  reached 
end  to  deposit  on  t«  inside  of  the  stapling  tube  or  on  the  Inner  veils  of 
the  isr i etc r,  rsther  thin  cr.  the  glass  slide;*  Shis  rail  lose  generally 
eaounta  to  less  then  3  par  cent  for  the  snail  droplets,  but  veil  loss  may  be 
as  tauch  es  10  or  even  £C  per  cent  for  drops  of  SO  microns'  diameter.  The 
extent  of  well  loss  should  be  determined  experimentally  for  the  specific 
test  conditions  imposed* 

Cor.cluslcn  P.ecardlng  The  Cascade  Ir-.pactor 

tten  dealing  with  liquid  dispersions  In  which  all  droplets  are 
under  50  xicrons  in  diameter  and  when  isokinetic  sampling  can  to  employed, 
the  cascade  lspactor  Is  an  excellent  instrument  for  deterndnlng  drop-size 
distribution*  The  suitability  of  the  instrument  decreases  as  maxi  mm  drop 
size  increases, owing  to  the  deposition  of  large  droplets  on  the  Inside  of 
the  sampling  tube. 


pyA.tics  c?  crjPisTS 

ace  un#-d -stance -velocity  re xauxons  cu  a  l— rup  depend  upon  two 
factors,  (I)  the  initial  velocity  of  the  drop,  ar.d  (2)  the  drag  of  the  sur- 
rouniing  air  into,  which  the  drop  is  injected.  Vor  large  drops  the  initial 
■velocity  dominates  the  dynamics  of  the  travel  for  an  appreciable  length  of 
time,  but  small  droplets  attain  the  »lr  velocity  in  an  extremely  short  interval 
of  time* 


-3-9- 


T.'rten  a  droplet  moves  relative  to  the  surrounding  air,  a  resistance 
or  drag  causes  the  droplet  to  accelerate  or  decelerate  until  its  velocity 
approaches  that  of  the  air.  iihen  the  particle  i3  at  rest,  v-ith  respect  to 
the  air,  the  drag  force  is  aero,  but  resistance  is  always  present  whenever 
the  particle  is  in  nation  with  respect  to  the  air.  This  r.ay  he  express e'd 
tv  “the  basic  equation 

R  *  X(v,v) 

where  R  is  the  resistance,  v  i3  the  difference  in  velocity  between  the  parti- 
el «  and  t,h«  medium,  and  v  is  the  acceleration  of  the  drop. 

The  exact  evaluation  of  this  function  is  complicated  in  the  case 
of  the  motion  of  droplets  ejected  at  high  velocity  from  an  atom! ting  nettle. 
The  functional  form  for  crag  is  complex  and  changes  markedly  over  most  of  tho 
range  of  interest  for  typical  sprays j  however,  in  the  limited  range  where 
the  droplets  have  slowed  down  and  are  moving  7.1th  laminar  motion,  the  simple 
Stokes'  law  applies.  Consequently,  only  brief  reference  v.ilj.  be  made  to 
Stokes'  law,  vdth  special  emphasis  on  its  limitations,  after  which  the  more 
complicated  functions  for  drag  will  bn  reviewed* 

Stokes '  Law 

• 

The  basic  St  okas'  lav;  equation  is; 

E  ••  onprv 

where  p.  i3  the  viscosity  of  the  fluid,  r  is  the  radiuo  of  the  r.  urical 
particle,  and  v  is  the  velocity  of  the  paiticle  relative  to  the  dium,  all 
in  cgs  units. 

One  of  the  principal  applications  of  Stokes  '  law  Is  the  evaluation 
of  the  terminal  velocity  of  freely  falling  p articles.  A  constant  terminal 
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velot  is  achieved  when  the  fluid  resistance  becorr.es  equal  in  magnitude 
and  t  rssite  in  direction  to  the  resultant  of  all  the  other  forces  acting 
on  the  particle.  Under  the  force  of  gravity,  the  equation  for  terminal  veloc¬ 
ity  is’>l*t 

2  (p-p‘)  r2  g 

e.  — — ■— «— -* 

.  *  9  n 

where  p  and  p»  are  the  densities  of  the  droplet  and  the  medium,  respectively, 
and  g  i3  the  acceleration  of  gravity, 

Stokes’  law  applies  accurately  only  for  particle*  moving  with 
laminar  motion  which  requires  that  the  Reynolds  number  for  the  particle  lie 
between  10“**  and  1,  Ae  a  consequence,  the  practical  range  cf  application 
for  droplets  of  unit  density  is  about  1  to  60  microns^.  Above  this  size 
special  empirical  equations  have  been  developed  to  determine  the  terminal 
velocity  of  droplets  up  to  1000  microns**. 

Figure  It  shows  the  terminal  velocity  of  water  droplets  for  diameters 
of  3  to  1C00  microns,  which  is  the  size  range  of  greatest  interest  in  the 
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Table  1  shows  v&luea  of  the  terminal  velocity  and  of  the  11  half-time » , 
(tj,),  for  five  sizes  of  droplets.  ’’Half-time"1  is  the  tine  required  for  a 
particle,  starting  from  rest,  to  attain  ons  -half  of  the  terminal  velocity. 

For  practical  purposes  small  droplets,  under  20  microns' diameter, 
attain  terminal  velocity  as  soon  as  they  begin  to  fall  freely  in  air,  and 
the  distance  traveled  before  attaining  terminal  velocity  i3  insignificant. 
Consequently,  the  dynamic  behavior  of  such  droplets  is  v/holly  dependent  upon 
the  movement  of  the- air  stream  in  which  they  fall,  J.arge  drops  fall  several 
feet  before  the  constant  terminal  velocity  is  closely  approached. 


Aa  mentioned  previously,  the  droplets  of  a  typical  spray  are  seldom 


raving  ritlt  laminar  notion,  as  required  by  Stokes '  law,  except  near  the  end 
of  their  trajectories*  Consequently,  the  taore  complex  and  practical  case  of 
droplets  roving  at  Reynolds  numbers  greater  than  1  will  be  considered,  priii  * 
_ clpalTy  from  t.hi»  standpoint  of  drag  ieient.  _  _  _ 


TABLE  1.  TERI-lL'PAL  VELOCITIES  AMD  HALT-TIME 
FOR  WATER  DROPLETS  IN  AIR 


T>ronl*t 

.tilouuj 

Terminal  velocity 

t ••  \  om  /•*« 

• 

Half-time, 

(t.  V  nee 
»  !*•  • 

20 

i.i 

0.0003 

60 

18 

O.cm 

200 

30 

O.cw 

500 

2?0 

o.i5 

1000 

500 

0.23 

Drag  Coefficients  Of  Drops 

The  drag  coefficient,  Cjj,  is  defined T  as; 

Cn  -  Srag 

(Frontal  Area) 

nhei6  gc  is  a3  defined  in  the  reference. 

For  steady  motion  of  droplets,  Cq  bscoir.es  a  function  of  Ro  alone, 

what* 


d  being  the  diameter  of  a  sphere  of  volume  equal  to  the  drop,  and  p'  and  p.' 

being  the  density  and  viscosity  of  the  air. 

For  droplets  Roving  with  streamline  Stokesian  motion,  the  expression 

for  drag  coefficient  is  extremely  simple j 
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Crs  ■  ~  (streamline  notion), 

***  Re 

However,  for  intermediate  nation,  and  for  turbulent  I.'ev.'tcnian  notion,  the 
expression  for  dreg  is  a  more- complex  function  of  Reynolds  number, 

Kughw  and  Gilliland7  recently  published  an  excellent  review  of  the 
motion  of  drops  in  which  they  awyiarized  the  available  data  on  drag  coeffi¬ 


cients. 

Figure  5  3 hows  the  drag  coefficients  for  spheres  and  drops  as  a 
function  of  Re  from  10  to  10,000, 

Kiighf^  and  rdlliland  point  out  that  the  analysis  of  the  notion  is 
greatly  complicated  by  the  period  of  acceleration  and  by  the  fluid  nature 
o i'  the  drop.  The  two  conventional  methods  for  determining  drag  coefficient 
are  (1)  direct  determination  of  drag  in  a  wind  tunnel,  and  (2)  the  measure- 
men^cf  terminal  velocities  of  falling  todies.  Because  such  measurements 
are  made  under  conditions  of  steady  flow,  the  important  effect  of  acoclera- 
ticr.al  drig,  which  generally  exists,  is  r.ot  included. 

Distortion  of  Drops 

As  shewn  in  Figure  5*  the  drag  coefficient  for  fluid  drops  is  not 
the  same  as  the  drag  for  solid  spheres,  at  the  higher  velocities,  owing  to 
distortion  of  the  drop.  Distortions  are  of  two  basic  types;  (1)  those  of 
an  equilibrium  nature,  and  (2)  those  of  an  oscillating  nature  resulting  from 
vibrations  about  the  equilibrium  position.  Distortion  occurs  because  the 
pressure  exerted  by  the  surrounding  fluid  on  a  moving  sphere  is  net  uniform. 
However,  the  theoretical  calculation  of  the  shape  of  the  drop  is  r.ot  possible 
because  the  pressure  on  the  surface,  which  In  turn  depends  on  the  clnpe  of 

the  drop,  is  not  known  exactly,  even  lor  a  sphere. 

-hi- 
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FIGURE  S.  DRAG  COEFFICIENTS  OF  SPHERES  AND  DROPS 
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(Hughes  and  Gilliland) 


Initial  Velocity  Of  Droplets  leaving 
CcntriFugi!'  Spray  ^osslcs 

A  factor  of  practical  interest  in  the  dynamics  of  sprays  is  tha 
initial  velocity  of  the  droplets  as  they  leave  tne  nozzl”.  Fcvikov^  r.sJe  a 
atiuiy  of  Luc  vciucity  of  liquid  ua  it  passes  threuga  the  atari 
charfcor  of  a  centrifugal  a  tori  zing  nozzle.  By  a  ccn3ider«tion  of  Bernoulli's 
aquation  ana  of  the  continuity  equation,  Fovikov  wss  able  to  develop  ex¬ 
pressions  for  the  velocity  and  the  thickness  of  the  film  of  liquid  issuing 
from  the  centrifugnl-type  nozzle. 

Ciffen  and  liuras  ze« 10  made  a  stroboscopic  study  of  the  discharge 
velocity  of  a  centrifugal  nozzle  ritfc  an  orifice  diameter  of  O.Oh29  inch. 

Tfc?  initial  discharge  velocity  of  the  fuel  spray  increased  from  U3  ft  per  sec 
at  ICO  psi  to  68  ft  per  sec  at  3CO  psl.  The  penetration  distance  v.-as  about 
27  inches  at  both  pressures. 


Dynamic  Effect  Of  Droplets  Dr.cn  Each  other 


Tinaliy,  mention  should  be  made  of  tlie  fact  that  the  dynamic  be- 

f  letj  ^  nf  ^11  niimo4  V>y  +  la  ^Tyni  m!  T,r  «d  ^  'Ir/.i'lnfo 

If  the  void  fraction  of  the  spray  is  85  per  cent  or  less,  the  stokes*  velocity 
13  reduced  by  more  than  $ 0  per  cent11.  However,  this  effect  of  other  r.oniby 
droplets  is  negligible  when  the  per  cent  of  voids  i3  greater  than  9?,  "»j  is 
true  a  short  distance  from  the  nozzle  orifice. 

Therefore,  in  a  fuel  spray,  the  effect  of  surrounding  droplets 
upon  the  behavior  of  any  single  droplet  need  be  considered  only  as  the  liquid 
first  emerges  from  the  atomizer.  1  fev;  centimeters  from  the  nozzle,  the 
concentration  of  droplets  is  so  dilute  that  each  droplet  r.v»  ie  cor.  v  id  ere  i 
to  tehavn  independently  of  th-  other  dro'letu  in  the  spray. 


Sm-— try  Of  Dynamics  0?  Proplets 


'f  he  rovenent  of  an  individual  droplet  d  3  3 uir.g  f To n  a  spray  nooolc 
dcpanas  upon  tba  initial  v*loclty  or  tha  drop  and  the  a  rag  aue  to  tb*  sur¬ 
rounding  air.  Drag  coefficients  of  reasonable  accuracy  *r«  available  in  the 
literature!  fco*ever,  the  an.Mysij  of  the  notion  is  greatly  complicated  by 
the  period  of  acceleration  and  by  the  distortion  of  the  drop* 

Additional  experimental  data,  especially  Hgh-speed  photography 
of  urcple la  ejected  it. ‘u.gh  veluwl v1(m> j  aic  i‘c^uli^<c<l  r^i  <i  uawvcx 
of  the  dynamic  behavior  of  droplets.  The  spinning-disk  atoniier  is  suggested 
as  an  excellent  tool  for  use  in  studying  the  trajectory  and  velocity  of 
droplets  • 
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L3QP-SI22  BiamiBU'flOUS 
hbl>  TU21R  AVsSi JWaS 


by 

Ralph  fa.  Thcaae 

V.£IGH~-.P  DISTRIBUTIONS  IKE  ’./SIOHTED  AVERAGES 

In  the  current  literature  there  are  several  types  of  distributions 
uhich  are  of  importance  In  describing  the  size  consist  of  sprays.  These  dis¬ 
tributions  arise  naturally  in  the  course  of  experimental  investigations  em¬ 
ploying  different  techniques.  For  example,  microscopic  exaninaticn  usually 
leads  to  droplet-size  distributions  by  frequency  ca-  number;  sedimentation  and 
impaction  lead  to  volume  (cr  weight)  distributions;  and  light  eatinction  leads 
to  surface  distributiens.  In.  actual  practice  the  true  surface  distribution 
and  true  volume  distribution  are  seldom  employed.  Instead,  the  surface- 
weighed  size  distribuUon  and  volume-weighed  size  distributiens  are  nor  a 
frequently  encountered,  hith  the  weighted  distributions  the  a -amis  is 
always  measured  in  terns  of  droplet  size  so  that  one  eaa  transform  fren  sur¬ 
face-weighted  distrituticn  to  volume-weighted  distribution,  for  uxanplo,  with¬ 
out  transforming  the  x-axis  iron  surface  to  volume  units.  Mathematically,  the 
surface-ueighted  size  distribution  nay  be  '.Tittcn  in  terns  of  the  surface  in- 
■crciuditis  mgr  ~ ~ -~i  ,  Thus*  ii'  rtA ^ v j ilojictss  tiio  o£ 

droplets  in  the  range  x  to  x  +  dr,  then  the  correci.tnding  increment  of 
surface  is  given  by 

x2  in(x)  dx 

dS(x)  =  — -  - 

j  X*  :n(x)  dx 
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cri  cumulative  distribution  function  of  the  sir  f  ac  e  -ue igh to d  sizo  distri- 
tutich  is  given  by 

S(x)  -  f"  d3(t) 

ii.-iilar  expressions  nay  bo  obtained  for  the  volume-weighted  size  distribution. 

xshic  i  gives  tic  curhiitivs  distribution  functions  for  the  various 
weighted  end  unweighted  distributions  usually  encountered  in  droplet  statis¬ 
tics.  Both  the  continuous  and  discrete  forms  are  given. 

Host  of  the  averages  ana  mean  droplet  diameter s  currently  employed 
in  droplet  statistics  are  readily  obtained  by  computing  tho  moments  of  the 
various  weighted  size  distributions..  For  example,  the  first  moment  of  the 
sixfaco'tcighted  size  distribution  is  given  by  the  expression, 

f  t(t*  fp{t)dt) / J  t2  fu(t)dt 

-o  /  -00 

cv  in  discrete  teras  by 

i  .  . 

£:  *i2  fn(xl)AXi 

1»1 


and  this  value  is  often  called  the  Saut-er,  or  "volurra-sui  face"  r.ean-droplet 


.j  _ _ _ 


^  .J  r,«  AA 


xributicn  may  bo  computed.  This  value  is  usually  called  the  weight-weighted 
mean  diameter.  From  the  calhematical  expressions  it  is  clear  that  the  various 
moments  of  the  weighted  size  distribution  are  expressible  in  terms  of  the 


»*>•>»•  -»,n  wwotavi 
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My  bo  vrltt»n  os  the  fi.-st  nmont  of  tha  tvrfaco-vol’.httfd  site  distribution 
er  as  tho  ratio  botecnn  tho  third  and  aeecnd  ntr.onto  of  tho  -inveichtod  olto 


tt  fallwi  that  tho  tarlcus  ntncnts  of  ih*  iIm  ilUtriUtUn  *f»  non* 
decreasing  as  tho  order  of  tho  nesent  inoroanea.  Thus,  tho  arlthnetio  man 
of  the  un'.fotfchtod  *i>a  distribution  la  lean  than,  cr  equal  to,  tho  ArltUuiitLe 
nean  of  tho  lcncth-volghtod  olio  dlatr Ibuticr.,  ‘.(hUh  in  turn  is  loos  than,  or 
equal  to,  tho  nrlthr.otlc  noon  of  tho  surfaca-i.‘oi£,htcd  alro  distribution,  oto. 
Tha  ncra  tho  droplot  alrca  aro  dloporood,  tho  stronG^r  Inequalities  will 
bo  bolifeon  tho  various  nor.nnta,  For  thlt*  rasen,  tho  ratio  of  tscr*,o  hlghor- 


crGor  suamt  to  a  louor-crdr*  nc^ont  la  often  tiV.en  aa  a  isoaauro  of  tho  dia- 
paralen  of  tho  dreplot  elien.  For  oxar.plo,  tho  ratio  of  tho  t'eipht-volGhtod 
noan  drop  alto  to  tho  volur.o-  cur  face  ir.ean  drop  alto  nlnua  coo  la  called  tho 
ccofflciont  of  variation  of  the  turfaeo-vo.lchtsd  olio  distribution.  In  a 


uiT(ti>*r  r<i?r«op 


to  tho  ur.voiplitod,  lontth-vclGlitod,  ar.d  voluiwrcichU’d  dlstrlt utlcns. 


Iho  variance  of  each  of  tha  veichted  dlutrlbutirna  rry  also  to 


i>>.j rosso  <i  In  terr.s  ef  tha  various  i-.cmnla.  ?cr  oxar.plo,  tha  variance  of  tha 
.virfacc-vaichted  r.lte  distribution  la  equal  to  tho  ]  redact  of  tho  third  and 
scscr.d  p.c^r.ts  of  *hf  si.29  ‘iijlrib'iticf)  nlw?  th?  •;fj»*in,*»  nf  itm 

second  swiit  of  tho  ufr.*ci(;htvd  dtstr Itu lien. 

;  ary  of  th.ro  equal!  t.lnr.  have  proved  i'.r.crtant  in  c*;.'-rl~-.:ital  in* 
v.-..'.l'.ntlcna  ir.vc'.v'. r.-  „ ;rf..cu* volur.e  relations,  Cf  particular  lnU.ro  at  la 
•  C*c*  tlot.  r.f  tetflr  inversely  proper*  1  anal  to  first,  ann’iil  ef 


the  surface -weighted  size  distribution,  and  this  moment,  in  turn,  i3  equal 
to  the  harmonic  mean  of  the  volume-weighted  size  distribution.  The  fact 
that  the  volume-surface  mean  diameter  is  greater  than  the  arithmetic  mean 
of  the  unweighted  size  distribution  and  at  tho  same  tame  is  less  than  tho 
arithmetic  mean  of  the  volume -v.’eighted  distribution  accounts  for  much  of 
the  acceptance  of  this  average  in  droplet  statistics. 

0‘hor  averages  of  importance  arise  from  the  following  typo  of 
consideration.  Calculate  the  diameter  of  a  hypothetical  droplet  whoso 
surface  area  is  equal  to  the  total  surfaco  area  of  tho  spray  divided  by  the 
total  number  of  droplets  in  Die  spray.  This  i3  another  kind  of  surface 
average  and  is  easily  shown  to  be  equal  to  the  go erne  trie  mean  of  the  mean 3 
of  the  length-weighted  and  unweighted  size  distributions.  Averages' of  tills 
typo  can  always  bo  reduced  to  a  geometric  mean  of  a  set  of  moments  of  the 
various  weighted  size  distributions. 

Table  2  summarizes  the  formulas  for  the  various  averages  usually 
encountered'  in  droplet  statistics. 


SPECIFIC  DISTOIHUTtOHG  KRh'utJSHTLY 
'EnPLOYr/iTlK  i'ROPlM’  STATISTICS 


Kcrrr.al  Distribution 


If  tho  diameters  of  the  droplets  in  .a  spray  have  a  probability 
density  function  given  by 


yn(x>  =  on  VZnT"  e  <  x  <  ®  , 

then  tha  droplet  sizes  are  said  to  be  normally  distributed.  Although  the 
normal  distribution  is  extremely  basic  to  most  phases  of  statistical  .analysis, 
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T/.BLF  2.  FOivrJLAC-  OF  ILlFCr.Wi.CE  I.*:  DROPLET  STATISTICS 

I  ARITK.IETIC  MEAX3 

Arithmetic  Mean’  of  Unw sighted  Size  Distribute on 

CD 

J  x  fn(x)dx 

For  continuous  distribution:  n  =  _ _ _ 

J  fn(r)dx 


For  discrete  distribution:  x. 


]>.  xifntxi)&xi 

i  =  1 _ 

fn(xi)ftxi 

1  =  1 


Ml  is  the  first  moment  of  the  continuous  size  distribution. 
xn  is  the  first  moment  of  the  discrete  size  distribution. 
xn  is  called  the  mean  dlomoter. 


Arithmetic  Mean,  of 'Length-', Velr.hted  Size  Distribution 


J  x  (*r„(x)dx) 

For  continuous  distribution:  u*\.  =  — _ 

*  I  CD 

f  xf„(x)dx 


For  discrete  distribution:  x. 


2  *1 (x(fn(x<)AxJ 
1  =  1 

CD  —  . 

*E.'  xifn(xi)Axi 
i  =1 


Min  is  the  first  moment  of  the  continuous  loncth-welghted 
_  size  distribution. 

*1  l“  Lbe  first  moment  or  the  discrete  lencth-weichted  size 
_  distribution. 

Xj  is  called  the  length-weirhted  mean  diameter.  . 


TAFLE  £.  FORMULAS  OF  IMRCRTAKCE  IK  1/iiOFLET  STATISTICS  (Continued) 


A_r  1 1  hr^e  tie  I  .•■■an  of  Surfaee-y,;eit’h:e-i  Size  Distribution 

CD 

/  x  (x2f,.  (z)dx) 

t 

Tor  continuous  distribution:  /u^,3  - 

J  x2fn(x)dx 

-<T) 

CD 

*i  I’' t2 rn (xi) '^i) 

Tor  discrete  distribution:  x3  =  -i— i - 

2.  *i2fn(xi)Axi 
i  =  1 

i*.’.  in  ♦.>»  fimt  moment,  of  the  continuous  ourfncc-wc-irht  ed 
•  x'  o 

size  distribution. 

x3  is  the  First  moment  of  the  discrete  surfacc-vioichted 
size  distribution. 

xg  is  called  the  Sauter  or  volume-surface  mean  diameter. 

Arithre ti c  ?‘cin  _o f_Vo ilurw-_{o r  y.'ol<th_t-)  '.'.'eirhted  0 i ze  Distri¬ 
bution 

/  x(x:'fn(x)dx) 

I  -CD 

For  continuous  distribution:  Mi»v  =  — - — - 

j  *®fn(x)dx 

-CD 

2.  5li(xi3 fn(xi)Axi) 
For  discrete  distribution:  xv  =  ---T-l - 

*  CD 

2..  Xi2fn(x)Axi 

i  =  1 

v  is  the  first  rr.cr.cnt  of  a  continuous  vol  un.e-v/oi^hted 
size  distribution. 

xy  is  the  first  moment  of  o  discrete  volume-weighted 
_  size  distribution. 

xy  is  colled  the  we i^ht-we ichted  mean  diameter. 

II  BASIC  II&fUAtrriFS  AMf'!!0  Th’^VAniOh".:  ARniJliTIC  tf.-ASS 
H.rrlovyp  l)f  Dh'or-'h'if  STATl's'flcY' 

For  continuous  distributions:  n  £  a  -  h\f  v 

For  discrete  distributions:  xn  <  xj  <  xs  <  xv 
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For  dlncreto  distribution: 


TA3I.Z  2.  FCnKULAS  CF  Ii.i?C3TAI.:CE  IT  Cr.OPLIT  STATISTICS  (Continued) 


IV 


Volui'.e-V.’eirhted  Size  Distribution 


For  continuous  distribution:  aS-  ~  — 


J  (x  -  /l.  v)2x"fr,(x)dx 


/  x3fn{xMx 


“00 

05 


For  discrete  distribution: 


2.  (xj  -  xv)2xi’fn(x1)Axj 
2  =  i  -  ] _ 

v  a) 

21  x  i3fn  (xi )  Axi 
i  -  1 


FORMULAS  SflOV/IKO  KFLATTOf.'S  SF/UV-FM  VARIAIICFS  AT  ID  AH  I  'ITtf.  TF  TIC.  fTKAfTS 

"  CF  sTsif  ;ii rrriuiiijf [OKI . . 

For  continuous  distributions: 

^r.2  '  l)  0‘1*  n)  "  (mi'.J2 
al2  =  l^l’s)  l'l)  "  (Mi'»i)S 
7s2  =  *  W»o)2 

OTy2  =  y)  (/i^,  y)  *  (/ipy)2 

For  discrete  distributions: 

®n2  =  *l*n  *  (*n)£ 
r’I2  =  xsxl  "  (xl2 


2  =  7  y 

It  1  '•  u 


, —  .  n  *  t 


In  tills  equation  ^  denotes  the  second  r'or.eiit  of  the  continuous 


tr/O  4  r*W»  r,A  »U»  JI  J»..  mV..- 

••  w  » tl 1  * ^  V-  *  V4  4  .>  V  i  iUUI/Atllt  J  JlUO 


v 


w 

/  x2(x3fn(x)dr.) 

-a: 

<X  ’  •*-— 

/  r.3fn(z)dx 


**  In  this  couation  xu  denotes  the  second  r/.rf  r.*  of  ?' 


k  Mi'pt' 


TABLE  2. 


FORMULAS  OF  If.^OHTAL.'CE  IK  DROPLET  STATISTICS  (Continued) 


I 


TABLE  2.  FOKJJLAS  OF  IMPORTANCE  III  DROPLET  STATISTICS  (Continued) 

VT _ FORMULAS _? OF_  THE  OLOMETRIC  MFAMS  OF  'IT?  VARIOUS  SIZE  DI3T3IBU- 

TlOt.'S  EIJCOUI.TEHED  It:  DROTLET  STATISTICS 


For  unweighted  size  distribution: 


—  i 

x„ 


rr  xt 

i  =1 


fn(xi)Axi 


1/'!>Z  fn(xi)Axi 

'  i  ~  1 


For  length-weighted  size  distribution: 


Jpj-  X(  axifn(xi>Axi'j  aXifntx^Axi 


For  surface-weighted  size  distribution: 


f®  ax12fn(xi)Axil  V^L-axi2  TnfxjJAxi 

*•  ’U  “  ]  1=1 

For  volume-weiGhted  size  distribution: 

f®  axi^^ri(xi'^xi]  3/i^*1oxi:!frl!xi)Axi 

J 

*  In  these  formulas  a  denotes  a  shape  factor. 

VII  FORMULAS  FOR  THE  HARMONIC  MF/n’tS  OF  TEE  VARIOUS  CIZ.E  LUhTfttD'J- 
Tio:ro  racou.ouifEb ~i: ' droplet'  statiuitco 

For  unv;eightea  size  distribution: 

CP 

-eL  '‘n(xi)Axi 

-  «  _  1  =  1 _ 

xn  "  cd 

*, 
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TABLE  FORMULAS  OF  IMPORTANCE  IK  DROPLET  STATISTICS  (Continued) 


VII  FORMULAS  FOR  TF.E  HARMONIC  MEANS  OF  THE  VARIOUS  SIZE  DISTRI3U 
TIOUS  ENCOUNTERED  IN  DROPLET  STATISTICS 


For  length- we ighted  size  distribution: 


x{'  =  1  =  1 


w 

21  xi^n  (xi)^xi 


-  xn 


2-  fn<xi)Axi 


For  6urfoce-v;eighted  size  distribution: 


-  »»■  -  i  -  1 
xfl  =  — 

6  co 


2  xi2fn(xi)*xi 


2T  xifn(xi)Axi 


For  volume-weighted  size  distribution: 


co 

y  x,"'fn(xi)Axi 
- _  i  =i _ .  - 

XV  -  00  ~  XB 

27  xisfn(x)&xi 

VIIT  RELATIONS  AMONG  AR1TC.IETIC ,  GEOMETRIC, AND_f!ASf.10MICJ.KAI^S_pF 

vARiou~s‘~3i'za  d i stImbTit ions 


General  Relat ion 

Harmonic  Moon  <  Goomotric  Mean  <  Arithmetic  Mean 
Unv;oightcd  size  distribution:  xn"  <  xn'  <  xn 
Length-weighted  size  distribution:  xj"<  x^'  <  Xj 
Surface-weighted  size  distribution:  x3"<  xa' <  xB 
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TA3LE  2.  FORMULAS  OF  IMPORTANCE  II!  DROPLET  STATISTICS  (Continued) 


VIII  RELATIONS  AMOI.'O ARITHMETIC,.  OEOISTRIC,  AND  HARMONIC  MgAKS  OF 
VARIOUS  SIZE  DISTRIBUTIONS 


General  Relation 


Volume-v.'eichted  size  distribution:  xv"  <  x^  < 
Mixed  inequality: 

—  1*.^  —  ,  1  —  —  M  —  |  —  —  ft  —  t  .. 

xn  1  xn  i  xn  s  X1  i  X1  1  X1  =  xs  i  xs  1  xs 

xs  is  called  the  volume- surface  mean  diameter. 
xy  le  caned  the  v;oi£ht-weichtea  mean  diameter. 


l.V  GEOMETRIC  AVERAGES  ENCOUNTERED  IN  DROPLET  STATISTICS 


The  diameter  of  the  droplet  whose  surface  ia  equal  to  the 
total  surface  of  the  6pray  divided  by  the  total  number  of  droplets  in 
the  spray  is  eiven  by 

Xgt  =  \fx^  %l 


The  diameter  of  the  droplet  whose  volume  is  equal  to  the 
total  volume  of  tho  spray  divided  by  the  total  number  of  droplets  in 
the  cprey  is  gives  by 


it  has  not  proved  satisfactory  in  the  representation  of  droplet  sizes, 
except  in  instances  vrr.ere  the  droplets  are  produced  by  condensation,  pre¬ 
cipitation,  or  by  chemical  processes.  Theoretically,  the  droplet  size3 
would  bo  expected  to  tend  to  be  uoriually  distributed,  provided  the  re¬ 
sultant  drop  size  is  the  cumulative  result  of  a  very  large  number  of  very 
small  chance  influences  where  the  influences  act  independently  of  each 
other  and  each  contributes  the  same'increment  to  the  resulting  drop  size. 
Since  most  observed  droplet  diebributions  are  ncn-normal,  it  may  be  con¬ 
cluded  that  the  underlying  hypotheses  are  violated. 

The  fact  that  the  normal  distribution  i3  fundamental  to  statistical 
inference  v:hilo  tho  normal  distribution  is  relatively  unimportant  in  droplet 
statistics  may  suggest  that  statistical  inference  may  not  be  valid  in  the 
realm  of  droplet  statistics.  However,  this  is  definitely  not  true.  If  an 
investigator  finds  that  his  data  are  normally  distributed,  he  may  apply  tho 
methods  of  statistical  inference  directly;  if  he  finds  that  his  data  are 
ncn-normal,  then  the  data  may  usually  bo  transformed  into  a  normal  form. 
Although  t.hn  unviflight.fid  droplet,  si  7.0  di  st.rihnt.i  on  is  non-normal  ;  it,  mnv 
happen  that  the  surface-weighted,  or  volume-weighted  distribution  may  be 
regarded  as  normal. 

Table  3  shows  tho  form  of  various  normal  size  distributions  en¬ 
countered  in  droplet  statistics. 

Log-i'ormal  Distribution 

In  many  instances,  a  ncn-normal  droolet-sise  distribution  can  bo 
transformed  into  a  normal  distribution  by,- replacing  tho  droplet  di ru-.o tor 
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TABLE  3.  VARIOUS  lOFJ.'.S  OF  THE  KOTJL5AL  SIZE  DISTRIBUTION  ENCOUNTERED 

IN  DROPLET  STATISTICS 


uith  thcbfsrltlw  of  the  drcplet  disraUr.  lit  this  cih  ths  droplet  sise 


Is  s*14  to  folic**  •  lc6*nor*»l  distribution  given  by 

„  (lcf  x  -  In*  i^)2 


x  v  r»  x 


“T* 
*lji  X 


x  >  0 


whore  lee  ia'  device  lecarltta  of  ttao  arlthratle  m*a  of  loc  x*  end 
"iofix  denotes  the  v4rl.1r.eo  of  lef  x. 

Theoretically.  tho  log-r.orrjl  diatrlbutlcn  can  bo  deducod  by  48* 
tURlnc  that  droplet  size  la  the  result  of  a  largo  r.urbor  of  nsall|  lndo- 
pendant  ir-xiltes  vhero  tho  t.  feet  of  each  Upelso  lo  proportional  to  tho 
alto  of  tho  dr  col*  t.  On  thlo  baole  alono,  It  soon*  clear  that  tho  loj- 
norrul  distribution  would  find  uidjr  applicability  In  dreplst  statistics 
than  tho  nomal  distribution. 

Table  It  chcua  tho  fora  of  tho  various  l0£-norrtal  distribution* 
eneamtored  In  drtplot  statistics. 


aoaln-Bnralcr  Distribution 

tho  Rcain-Rsanlor  diatrlbutlcn  la  a  wclcht-welghted  size  distrl¬ 
butlcn  that  has  been  successfully  srpllod  to  sir*  ccn:lcta  reizltlng'  frea 
(rindlnp.  !latho;iatl cully,  the  probability  density  function  Is  £lvon  by 

y„{x)  *  nhx“"1  •*h9,a  f  x  >0  , 

*.: hc.-a  n  and  b  are  r.ur.erleal  paranctcra.  Tha  parnrotrr  n  la  a  characteris¬ 
tic  of  tbs  substance  of  which  the  particles  are  compared;  tho  vaJ  jo  of  b 
increor-s  »*♦■*»  j«rro»*in*  fineness,  the  '•ettht*’"*!  ««-»*»  site 
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IV.  Yc  1  ur.e -v.-e  i  g  ht  ed  log  normal  distribution: 


(lcr  x  -  log  xvV 


Vy  I*)  “  , 


r*  (7n  2 

log  X  ,  X  >0 


whero 


log  xv'  =  J  {log  x)  nx^lxjdx 
Jo 

*> 

^log  X2  s  J  *  log  x  -  log  xv’)2  ax3fn(x)dx  * 


ggtgagsgas? xssvsjy a errgnn  jxtx  sagf- aur aagggaaaaBBr .  •  t»rjr- 

*  In  these  formulas  a  denotes  an  appropriate  shape  factor. 
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is  easily  found  to  be  equal  to  l/>Vb  r(l/n  ♦  1),  where  r(l/n  ♦  1)  re¬ 
presents  the  generalized  factorial ,  Frcr.  this  it  is  clear  that  the  greater 
the  fineness  of  the  .consist,' the  smaller  the  weight-weighted  mean  diameter. 
The  Rosln-Rairmler  distribution  can  be  often  fitted  to  data  that  are  too 
skewed  to  be  fitted  by  a  log-normal  distribution. 

Table  5  shows  the  various  forms  'of  the  Rosin-Rammler  distribution. 


Muklyana-Tanasawa  Distribution 

In  correlating  the  data  obtained  from  air  atomization  of  liquids, 
Hukiyama  and  Tanasawa  used  the  following  unweighted  size  distribution! 

yn(x)  "  |y  x2  e-*>*  ,  x>  o  . 

This  equation  is  one  member  of  a  family  of  statistical  functions  known  as 
the  Gamma  distributions. 

Table  6  shows  the  various  forms  of  the  Hukiyama- Tanasawa  Distri¬ 
bution. 


Other  Dis  tributions 


In  many  instances  the  preceding  size  distributions  will  not  yield 
an  adequate  fit  to  the  observed  data  so  that  other  distributions  must  bu  con¬ 
sidered.  Important  alternative  distributions  include  the  fcllcwingi 
(l)  lieibuli's  Vclume-v.eighted  Size  Distribution  is  given  by 

vn*>l 


yn<*>  -  = 


-  /*  "  x,\n” 

?v\  3cv  / 


6  N  xv  ■/ 
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TABLE  5.  VARIOUS  FORMS  OF  TEF  ROSIK-RAJ.1XER  DISTRIBUTION 


1.  V.'eieht-v/eighted  size  distribution: 

yw(x)  =  n  bxn_1  e_b  x°  ,  x  >  o 

Uelcht-weighted  mean  droplet  size  =  b"*/n  Tfl  ♦  1/n) 

4l 

II.  Surfooe-v.eighted  si2e  distribution 

yB(x)  =  a  nb  xn"  2  e  _b  x"  ,  x  >  0 
Suriaco-we tinted  mown  droplet  ei»e  -  a 

III.  Unweichtcd  eizo  distribution  * 

yn(x)  =  anbxn"  *  e  "b  x°  ,  x  >  o 

Mean  droi  iet  diameter  =  ab^/n  P  (X  -  8/h) 

*  In  these  equations  a  denotes  s  shape  factor. 


TABLE  6.  VARIOUS  FORMS  OF  THE  NUKIYA!.!A-TAJ;A$A’..'A  L[1  ST?.  I  BUT  I  ON 


I.  Unweighted  size  distribution: 

yn(x)  =  ~  x2  e'bx  ,  x  > 
2! 


Unweighted  mean  diameter:  3/b 
II.  Surface-weighted  size  distribution:' 

,  X  >  o 


ab3 


y°M  '  IT  x  8 


Surface-weichted  mean  diameter:  5ta/?’j3 


III.  Volume-weighted  size  distribution: 

yvl*)  s  ^rr  x®  *  "b*  ,  x  >  o 


21 


Volume-weighted  mean  diameter:  6!a/2b4 


*  In  these  formulas  a  denotes  a  chape  factor. 
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where  x0  is  the  smallest  reserved  droplet  size.  This  distribution  has 
properties  sir.ilar  to  these  of  the  P.csin-F.ammler  distribution. 

(2)  The  Grar.-Oharlier  Series  represents  an  arbitrary  droplet- 
size  distribution  as  an  infinite  series  expansion  in  terns  of  the  Hermite 
polynomials  as  follows J 

yfx)  =  ak  , 

where 

00 

ak  =  k.'  J  Hn(x)  y(x)  dx  > 

m 


^(^(x)  denotes  the  k**1  derivative  of  tho  error  function,  and  Hn(x)  denotes 
the  Hermite  polynomial  of  decree  n.  For  data  which  are  approximately 
normal,  the  first  few  terms  of  the  scries  expansion  .are  usually  sufficient. 
(3)  The  Gamma  Distributions  are  given  by 

yW  "  r<«  .  i)  a®*1  xa  '■x/a  •  x  >  0  1 


•  M..M . .  >r-  J4  J  -  •  P  .  ..  .t.c  .  _  /•  it  «  /* - j  t ..  -  i* 

UiD  jumijai  la*  ianuoiuwa  ^iiuutiun  io  a  p-*i  oiuuia  i  MCiiiubj  \j  I  w**»3  idi'iiijf  Ui 


Gamma  distributions. 

(h)  Any  non-negative  function  which  has  a  finite  .integral  over 
the  range  considered  may  servo  as  a  size-distribution  function,  However,  to 
enablo  one  to  make  statistical  inferences  tho  parameters  of  tho  chosen  distri¬ 
bution  must  have  known  statistical  distributions.  A  distribution  which  gives 
an  excellent  fit  may  not  be  satisfactory  for  the  purposes  of  statistical  in¬ 
ference  if  tho  parameters  do  not  have  known  distributions.  For  this  reason, 
considerable  care  must  bo  exorcised  in  employing  tho  various  curve-fitting 
techniques  which  introduce  new  parameters  into  the  distribution. 
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THE  CHI-SQUARE  TEST  OF  SIGNIFICANCE 


Having  chosen  a  particular  type  of  distribution  function,  the  con¬ 
stant  parameters  may  be  computed  so  that  the  curve  gives  the  best  fit  to  the 
data.  Graphical  analysis  is  often  helpful  in  determining  the  order  of  magni¬ 
tude  of  the  parameters.  The  usual  curvo-fitting  methods,  such  as  least 
80,uares,  are  generally  sufficient  for  this  purpose,  However,  analytical 
methods  should  supplement  the  graphical  methods  in  the  final  determination 
of  the  values  of  the  parameters. 

aI**#.**  4Vm  Jl  A.w  .i!  ...  i.».  .  » . i  • 

•  VIM  puk*  w*  W'i  *UI  uawvt  vWUWAVIl  AWUVVAWIIj  UIW  J.M  *  V  O  WJL  “ 

gator  should  test  the  distribution  against  the  observed  data  to  determine 
vhothor  the  deviations  between  them  can  be  attributed  to  chance.  To  ac¬ 
complish  this  the  following  expression  is  computed: 


j(2  £  (expected  value  -  observed  value  )2 

expected  value 


This  gives  the 'observed  value  of  Chi-Square,  If  this  value  is  very  large, 
it  is  extremely  improbable  that  the  observed  droplets  originated  in  a  spray 
having  the  oise  distribution  being  tested.  In  order  to  make  the  test  quanti¬ 
tative,  a  level  of  significance  is  chosen,  say  a,  and  the  corresponding  X2a 
is  obtained  from  the  X2  tables.  Then,  if  the  observed  value  of  exceeds 
the  value  of  X2-  ,  the  null 

m  *  ■  m . .  »■  ^  •  •*  •*/  *-* •  —  .  vi"  ■* 

distribution  should  be  considered. 

J.t  may  bo  noted  that  sinco  the  curve  giving  tho  boat  fit  from  a 
p-.u*'  W.U.ty  point  of  view  is  that  curve  which  renders  tho  valuo  of  X2  a 

This  .prlnciplo  may  bs  employed  t,o  compute  t.'uo  nur.-ricrl  values  of 

inters  to  bo  used  in  tha  distribution. 
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Best  Available 


It  should  also  be  noted  that  the  theoretical  distribution  nay  fit 
the  data  too  well.  In  this  instance  the  observed  X^  has  a  value  so  small 
that  the  probability  of  obtaining  such  a  value  is  very  close  to  aero.  This 
result  often  occurs  when  too  many  parameters  have  boon  introduced  into  the 
distribution  function.  In  general,  the  more  parameters  one  has  at  his  dis¬ 
posal,  the  better  the  distribution  will  fit  the  data.  By  eliminating  a 
parameter,  an  overfitted  distribution  may  berendered  less  exact  and  a  more 
reasonable  value  of  may  be  obtained. 

If  the  observed  value  of  &  has  a  probability  of  occurrence  between 
0.05  and  0.95,  it  is  usually  agreed  that  the  distribution  being  considered 
will  serve  as  a  basis  for  correlating  the  droplet-size  data.  In  some  in¬ 
stances,  more  than  one  distribution  may  give  a  satisfactory  fit.  In  those 
instances,  the  distribution  with  the  most  reasonable  physical  interpretation 
should  be  given  preference. 


Inadequacy  Of  Relative  Error 


A  simple  calculation  will  chow  that  the  customary  use  of  relative 
error*  as  a  measure  of  the  difference  between  theoretical  and  ebserved  re¬ 
sults  is  not  adequate  when  probability  considerations  are  involved.  If  the 
probability  of  obtaining  a  droplet  in  a  specified  size  range  is  denoted  by 
£»  and  if  a  sample  of  size  n  is  drawn  from  the  population,  the  expected 
number  of  droplets  in  the  given  range  is  equal  to  np.  In  general,  the  qb-  ■' 
served  number  of  droplets  will  bo  different  from  the  expected  number.  The 
ebserved  number  may  be  written  as  np*,  where  p'  denotes  the  empirical 


.*  Relative  error  is  defined  as* 


4- J  -ml 
ii  UAPUdk 


Theoretical  value 
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probability  of  obtaining  a  droplet  in  the  tiven  range,  The  relative  error 
between  tho  observed  rnd  o;«cUd  results  is  given  by 


r«P*  -  nP  . 
np 


t 


And  this  expression  Is  Independent  of  the  site  of  sarspio.  The  contribution 
to  X2  because  of  the  different*  beUnon  the  ebsurvod  end  expected  reeulte 
ie  equal  to 

■fog----  "Eli  .  n(p»  -  p)» 

np  p  » 

end  thie  value  ia  directly  proportional  to  the  also  of  tho  sanple.  Thua, 
doubling  the  sanple  site  doublet  the  value  of  X2  whenever  the  relative 
error  remains  constant.  Ccnsequen  tly,  the  probaoility  of  obtaining  the 
earns  relative  error  in  a  doubled  site  of  sample  is  reduced.  This  la  in 
accord  with  the  intuitive  notion  that  In  larger  samples  the  deviations 
due  to  chance  shculd  "settle  dejsj",  ar.d  to  eppese  a  constant  relative  error* 


Llrdtatlcna  Of  statistical  Analysis 

Statistical  analysis  is  nost  powerful  vhen  dealing  uith  chance 
errors.  In  aany  instances,  iKwovor,  tho  errors  nay  fcu  SjTIi^ter^tdC ,  and 
systenatie  errors  do  not  hovo  statistical  pc*:ulatims.  If,  for  example, 
droplet  ccunts  ore  nade  consistently  too  high,  or  weight  netsurtuvnta  ore 
Uide  consistently  too  lc*.i,  the  resulting  bias  In  the  sa.-plo  cannot  be  at¬ 
tributed  to  chance  influences,  in  this  instance  statistical  lusts  of 


significance  become  meaningless,  and  the  investigator  may  argue  that  his 
problem  is  thereby  reduced  to  cne  of  simple  curve  fitting  where  no  regard 
need  be  given  to  probability  considerations.  However,  it  should  be  em¬ 
phasized  that  unless  the  systematic  errors  can  ba  determined  and  their 
effects  eliminated,  the  fundamental  problem  of  droplet  statistics  can  not 
be  solved.  That  is,  the  investigator  will  not  be  able  to  state  the  drop- 
size  distribution  of  the  spray  from  which  the  sample  was  draw.  If  he 
contends  that  systematic  errors  cannot  be  eliminated,  then  he  must  admit 
that  the  parent  drop-size  distribution  is  both  unknown  and  unknowable.  On 
the  other  hand,  if  the  systematic  errors  can  be  evaluated  and  eliminated, 
then  a  complete  statistical  analysis  should  produce  a  distribution  which 
describes  the  parent  droplet  population  with  a  specified  degree  of  reli¬ 
ability. 

RET tea 
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?e view  cf  mass  aj.d  e.'efsy  Taa::sfe*  feiveen  liquid  duplets 
a:;d  surrcui.dtks  c-^es* 

Dy  Charles  C  •  Graves  and  0.  V.  Bahr 

National  Advisory  Committee  fcr  Aeronautics 
Cleveland,  Ohio 


Formulation  of  design  criteria  for  high  performance  Jet  engine  com¬ 
bustors  or  for  industrial  equipment  such  as  spray  dryers  and  cooling 
towers  requires  accurate  knowledge  cf  liquid  vaporization  rates .  Condi¬ 
tions  fcr  evaporation  vary  widely  among  these  various  commercial  and 
military  applications.  As  an  example,  contrast  the  relatively  low  gas 

4  m  Ima  *  ««4  4-U  •*<■»  4  f  4  r.n  i«m  ^  a  A  AA 

•  "  -  •••  - -  - —  -r  - -  r  — 

second  in  the  ramjet.  Or,  contrast  the  relatively  low  temperatures  in 
the  cooling  tower  with  conditions  in  the  turbojet  combustor  where  fuel 
is  sprayed  into  regions  approaching  stoichiometric  flame  temperature. 
Accordingly,  as  shown  on  the  first  slide,  this  paper  is  concerned  with 
conditions  for  single  drops  ranging  from  evaporation  in  lew  temperature 
quiescent  surroundings  on  through  to  a  discussion  of  the  effects  of 
forced  convection  on  heat  and  mass  transfer  and,  finally,  to  the  con¬ 
sideration  of  the  unsteady  state  process  existing  during  the  initial 
period  of  evaporation.  Since  another  review  paper  will  be  presented  on 
spray  drying,  the  final  portion  of  this  paper  will  be  limited  to  but  a 
brief  discussion  of  the  evaporation  of  sprays. 

» 

First,  let  us  consider  the  case  of  Bteady  state  drop  evaporation  in 
quiescent  surroundings.  This  case  is  of  interest  since  it  offers  a  simple 
Illustration  of  the  principles  involved  in  evaporation  over  the  ran<?e  of 
conditions  covered  in  this  paper.  In  addition,  the  calculated  evaporation 
rateB  may  be  used  as  reasonable  first  estimates  of  evaporation  rates  for 
some  applications. 

The  assumed  physical  model  is  shown  on  the  next  slide.  Vapor,  at  the 
saturation  pressure  fcr  the  drop  surface  temperature,  diffuses  to  the  sur¬ 
rounding  atmosphere.  The  drop  surface  temperature  is  maintained  at  a  tem¬ 
perature  sufficiently  telcv  that  of  the  surroundings  to  provide  the  vapori¬ 
zation  heat  transfer  requirements.  We  see  thi_n  that  the  drop  evaporation 
rate  can  be  expressed  in  terms  of  either  mess  or  heat  transfer  equations. 

In  the  derivation  of  these  equations  it  Is  assumed  that  a  quasi - 
stationary  state  exists  ouch  that  the  vapor  concentration  and  temperature 
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gradients  around  the  drop  are,  at  every  instant,  those  corresponding  to 
equilibrium  values  for  the  existing  drop  size.  It  is  also  assumed  that 
the  temperature  in  the  drop  interior  is  constant  at  its  equilibrium 
value  For  the  condition  of  spherical  symmetry,  the  differential  heat 
and  mass  transfer  equations  for  a  spherical  surface  a  distance  r  from 
the  drop  center  are  given  by  equations  (1)  and  (2)  In  equation  (l) 
the  term  c-  -  is  the  enthalpy  change  of  the  fuel  vapor  between 
the  drop  eu  and  the  surface  considered,  L  is  the  latent  heat  of 

vaporization  at  the  surface  temperature,  and  dw/dt  is  the  drop  eva¬ 
poration  rate.  This  equation  simply  equates  the  heat  transferred  by 
mass  movement  of  the  diffusing  vapor  to  that  transferred  by  thermal 
conductivity.  Equation  (2)  is  the  application  of  the  Stefan  diffusion 
equation  (ref.  1)  to  drop  evaporation. 

For  evaporation  in  quiescent  surroundings  Lheee  equatioua  ore  inte¬ 
grated  between  the  drop  surface  and  infinity  to  obtain  equations  (3)  and 
(4).  For  these  integrations,  cD,  k,  and  D/T  were  taken  as  constants 
and  represent  mean  values  between  the  drop  surface  and  the  surroundingo 
Vapor  pressure  in  the  surrounding  atmosphere  was  assumed  negligible. 
Predicted  drop  evaporation  rates  would  be  obtained  from  simultaneous  solu¬ 
tion  of  these  heat  and  mass  transfer  equations  along  with  the  assumption 
that  the  vapor  pressure  at  the  drop  surface  corresponds  to  saturation 
pressure  for  the  drop  surface  temperature. 


For  evaporation  in  low  temperature  ourround. ngs  the  sensible  heat 
change  of  the  vapor  between  the  drop  surface  and  surroundings  is  small 
compared  to  the  latent  heat  of  vaporization  Cons*  ;uently,  equation  (3) 
reduces  to  equation  (S)  on  the  next  slide.  This  equation  is  similar  in 


form  to  that  of  simple  heat  transfer  to  a  sphere  In  quiescent  surround- 
1  for  cVapOi'etwlOu  Wuci’c  uli<S  Vd  per  pressure  u*  «he  drop 


surface  is  small  compared  to  the  total  pressure,  equation  (4)  reduces  to 
equation  (6),  the  Langmuir  .  -uation  (ref  2)-  It  la  seen  that  equations 
(5)  and  (6)  are  of  similar  form,  illustrating  the  widely  known  analogy 
between  heat  and  mass  transfer  On  the  other  hand,  as  was  seen  from 
equations  (3)  and  (4)  on  the  previous  slide,  the  analogy  breaks  down  for 


evaporation  in  high  temperature  surroundings  In  either  case,  however, 
drop  evaporation  is  proportional  to  drop  diameter. 


Fquation  (6)  has  been  verified  for  evaporation  of  small  drops  in 
room  air  by  the  work  of  Topley  and  Why t.luw -Gray  (ref  3);  Houghton 
(ref.  4);  and  Whytlaw-Gray  and  Patterson  (ref.  5) 


We  have  seen  that  equation  (6)  applies  only  when  the  surface  vnpor 
pressure  is  small  compared  with  the  total  pressure.  Consider  for  a  moment 
that  other  corrections  might  be  made  to  equation  (6)  If  we  divide  this 
equation  by  drop  surface  area,  we  see  that  the  evaporation  rate  per  unit 
surface  area  is  inversely  proportional  to  drop  dimeter.  In  addition, 
this  rate  is  inversely  proportional  to  total  pressure  as  determined  by  the 
pressure  dependency  of  the  diffusion  coefficient  Consequently,  as  we 
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approach  stiff iei*»r.ti/  small  ir 
rate  per  unit  area  as  cl  tamed 
theoretical  maximum  evaporeti  ■*. 
tained  from  gas  kinetics  '‘ref. 
tton  rate  for  a  drop  evapcratl 
equation  (’)  where  la  1/4 
and  a  Is  the  evaporation  or 


and  tort*!  pr.-saurcs,  the  evaporation 
.  frer.  equal  sen  tv)  could  exceed  the 
n.  rate  for  evaporation  In  a  vacua-?,  as  ch¬ 
ic)  .  lata  theoretical  maximum  e  vapor  a - 
ng  iBOtr.err.aily  In  a  vacuum  is  given  by 
tr.~  arithmetic  mean  -clccular  velocity 
arcomcdatlcn  coefficient. 


The  case  of  evaporation  of  very  ena/i  drops  wtis  treated  theoreti¬ 
cally  by  Fuchs  (ref.  7).  He  considered  diffusion  to  start,  not  at  the 
drop  surface,  hut  at  a  distance  approximately  one  mean  free  path  from 
the  drop  surf ace  The  mass  transfer  drop  evaporation  equation  ac¬ 
cording  to  this  picture  la  given  by  equation  (8)..  Here  A  is  the  dis¬ 
tance  between  the  drop  surface  and  the  surface  where  the  diffusion  pro¬ 
cess  is  assumed  to  start.  For  drop  disasters  large  compared  to  A, 
this  equation  reduces  to  the  Langmuir  equation.  For  very  small  drops, 


~  4>U* 


rate  approaches  that  of  a  drop  evaporating  in  a  vacuuia  as  given  by 
equation  (7).  Similarly,  for  the  larger  size  drops,  thiB  theoretical 
maximum  evaporation  rate  is  approached  as  we  go  to  sufficiently  low 
pressures  where  A  is  again  large  compared  to  the  drop  diameter.  The 
Fuchs  theory  was  confirmed  experimentally  by  the  work  of  Bradley  and 
his  co-vorkero  (refs.  8,  9,  ani  10). 


For  most  applications,  the  correction  according  to  the  Fuchs  theory 
is  mainly,  of  academic  interest.  The  effect  generally  becomes  noticeable 
at  atmospheric  pressures  cnly  when  we  reach  drop  sizes  in  the  order  of 
one  micron.  Similarly.,  ccrrection  for.  the  effect  of  surface  tension  on 
surface  saturation  vapor  pressure  generally  becomes  noticeable  only  for 
drops  smaller- than  one  micron.  Finally,  Lu.,hak  and  Langstroth  (ref.  11)  \ 

have  shown  that,  errors  resulting  fr-.n  the  assumption  of  quasi-otatlonary  \ 
state  are  negligible  under  moat  circumstances. 


Now  let  us  consider  the  case  ot  evaporation  In  high  temperature 
quiescent  surroundings.  There  has  beer,  little  experimental  work  done 
for  this,  region.  All  of  thiu  work  appears  to  huve  been  done  for  the 
cose  of  single  drops  burning  In  quicieent  .-.i  -.rapheres.  Since  drop 
burning  has  been  treated  theoretically  as  a  special  cae  of  drop  evapora¬ 
tion,  a  brief  discussion  of  the  work  is  warranted:.  For  the  high  tempera¬ 
ture  evaporation  case  of  burning  drops  we  use  the  heat  transfer  equation  v 
and  avoid  use  of  the  muss  transfer  equation  for-  fuel  vapor  by  assuming 
the  surface  temperature  to  bo  the  fuel  lolling,  point..  Since  we  nre 
dealing  with  large  temperature  differences,  the  error  involved  In  this 
assumption  is  minor.  The  next  slide  presents  the  model  for  the  burning 
drop  as  proposed  by  Spalding  (ref.  1':).  In.  this  picture  oxygen  diffuses 
Trqm.  the  surrounding  atmosphere  C  to  a  burning  cone  of  negligilile 
thickness  located  at  B,  while  exhjuyjt  products  diffuse  In  the  opposite 
direction.  Heat  balance  equations  are  set  up  for  regions  AB  and  BC. 

In  addition,  c  mss  transfer  equation  is  set  up  for  the  diffusion  of 
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oxygen  from  the  surrounding  atmosphere  to  the  turning  zone.  The  drop 
turning  or  evaporation  rat*  U  obtained  from  simultaneous  solution  of 
these  equations  This  burning  rate  is  given  by  equation  (1),  which 
represents  the  simplified  equation  where  thermal  conductivity  and 
specific  heats  are  considered  to  be  constant  Accordingly,  these  values 
are  taken  as  mean  values  for  the  appropriate  regions  of  integration. 
Equation  (2)  gives  the  simple  heat  transfer  equation  for  low  temperature 
evaporation.  From  comparison  of  these  two  equations,  we  can  see  that 
the  term  in  the  brackets  of  equation  (l)  represents  an  effective  tempera¬ 
ture  difference  for  evaporation.  It  is  noted  that  for  this  case  of  high 
temperature  evaporation,  the  sensible  heat  change  of  the  fuel  vapor  is 
large  compared  to  the  latent  heat  of  vaporization.  Consequently,  serious 
error  would  result  by  neglecting  the  enthalpy  of  the  fuel  vapor.  This 
is  illustrated  by  the  next  slide.  On  the  ordinate  is  plotted  the  ratio 
of  evaporation  rate  predicted  by  the  simple  heat  transfer  equation  to 
uimi.  predicted  by  the  equation  accounting  ror  vapor  entnaipy.  Curves  are 
presented  for  iaooctane  and  water  drops  evaporating  in  quiescent  surround¬ 
ings  at  various  temperatures.  It  is  seen  that  at  the  higher  temperatures, 
the  simple  heat  transfer  equation  would  predict  much  higher  evaporation 
rates  for  the  case  of  iaooctane  drops.  This  effect  1b  much  lees  pro¬ 
nounced  for  water  which  has  a  high  latent  heat  of  vaporization.  Effects 
of  diffusing  vapor  enthalpy  changes  on  predicted  evaporating  rateB  were 
treated  theoretically  by  Ackermann  (ref.  13)  and  by  Colburn  and  Drew 
(ref.  14). 

Spalding  (ref.  15)  conducted  some  experimental  work  on  a  burning 

sphere.  However,  the  sphere  sad  a  aiwut'wcr  of  li  inches  and  consequently 

& 

was  subject  to  considerable  free  convection  effects  which  were  difficult 
to  predict.  In  addition,  it  is  believed  that  possible  high  radiation 
heat  transfer  to  the  sphere  may  have  obscured  the  results. 

Godeave  (ref.  16)  determined  the  burning  rates  for  a  number  of  hydro¬ 
carbon  fuels  from  photographic  measurements  of  silhouettes  of  burning 
drops  suspended  in  quiescent  air.  Initial  drop  diameters  were  approxi¬ 
mately  1700  microns.  An  approximate  relationship  was  found  between  burn¬ 
ing  rate  and  the  enthalpy  change  between  liquid  fuel  at  the  Initial  drop 
temperature  and  vapor  fuel  at  the  boiling  point.  Turning  rate  ms 
directly  proportional  to  drop  diameter. 

Codsave  (ref.  17)  also  treated  the  turning  of  single  drops  theoreti¬ 
cally  by  ueing  the  heat  transfer  equation  for  r.vapuratlon  in  high  tempera¬ 
ture  surroundings .  Predicted  burning  rates  were  obtained  by  assuming  a 
temperature  for  the  burning  zone  surrounding  the  drop  and  by  substituting 
film  thickneeses  for  the  region  between  the  drop  surface  and  turning  zone 
as  determined  from  photographs .  Good  agreement  between  predicted  and.  ex¬ 
perimental  values  was  obtained  by  this  method. 
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Graves  (ref.  16)  used  the  deleave  pr* led:. re  to  determine  the  effect 
cf  oxygen  concentration  on  th-  turning  rate  cf  n ingle  iococtano  drops. 
Beth  the  absolute  burning  rates  and  changes  In  burning  rate  with  oxygen 
concentration  were  predicted  rather  closely  by  equations  baaed  on  the 
Spalding  theory. 

Even  for  the  relatively  saa.ll  drop  sizes  investigated  by  Gcdsave 
(ref.  IS)  and  Craves  (ref.  13),  the  free  convection  effects  resulting 
from  the  high  temperature  gradients  nay  have  been  appreciable.  In  their 
paper  on  drop  evaporation,.  Ranz  and  Marshall  (ref.  19)  suggested  a  cor¬ 
relation  for  free  convection  similar  in  form  to  the  forced  ccnvectlon 
correlation.  If  this  relation  is  used  in  conjunction  with  a  modified 
Groshof  number  (ref.  20)  for  the  burning  drops,  we  n ‘ght  expect  buoyancy 
forces  to  Increase  burning  rate  in  the  order  of  50  percent  over  that  for 
on  infinite  film  thickness.  In  »*•«  l,h<»  drop  diameter  dependency  of 

the  Graahof  number  drop  burning  rate  would  then  be  proportional  to  a 
power  of  diameter  greater  than  unity.  However,  data  for  the  burning  drops 
did  not  appear  tc  show  this  trend.  Since  the  drops  were  initially  at  a 
temperature  considerably  below  the  equilibrium  value,  a  possible  factor 
offsetting  the  effect  of  free  convection  would  be  the  unsteady  state 
heating  of  the  drop  interior.  This  process  will  be  considered  later  in 
the  paper. 

We  have  considered  drops  evaporating  In  quiescent  surroundings. 

For  negligible  free  convection  effects  the  fill?  thickness  was  considered 
to  be  infinite.  For  this  case  we  also  assumed  spherical  symmetry  of  the  • 
transfer  paths.  Now,  let  us  consider  drop  evaporation  under  forced  con¬ 
vection. 

With  forced  convection  we  no  longer  have  spherical  symmetry.  In 
fact  there  may  be  a  largo  chnnge  !r.  local  transfer  rate  over  the  drop 
surface.  An  example  of. this  change  is  shown  In  the  next  slide.  Here 
local  trunsfer  rates,  expressed  as  a  fraction  of  the  transfer  rate  at  the 
forward  stagnation  point,  Is  plotted  against  the  angle  from  the  forward 
stagnation  point.  The  data  w.-re  obtained  ry  Frooasl Ing  (ref.  21)  for  the 
sublimation  of  a  naphthalene  bead.  Boundary  layer  theory  jr olio to  that 
the  local  transfer  rate  should.  be  a  maximum  at  the  forward  stagnation 
point,  decrease  gradually  tc  a  minimum  near  the  Depuration  point  and  then 
rise  gradually  to  a  uew 'maximum  ut  the  ivur  stagnation  point.  This  type 
of  chnnge  ir  0-  **n  to  be  followed  In  the  figure. 

Ran.  and  Marshall  (ref  ■  19)  determln.-d  the  temperatur  e  pr»  f  lie  around 
an  evaporating  water  drop  ns  shewn  in  th»  next  olid-.  It  la  seen  tiiat, 
for  a  Reynolds  number  of  approximately  R,  the  boundary  layer  thickness 
and  drrp  diiuncter  ore  approximately  the  nn’i:e„  1 'ivi  tiie  ehunge  in  spacing 
of  the  lines  r.f  constant  temperature  d i ff.-r«'iic«  it  la  also  aeon  that  there 
again  a  gradual  decrease?  in  local  trarjsfor  rale  from  the*  value  ’ 

at  the  forward  stagnation  point. 
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Correlations  if  experimental  data  neglect  char.g-.-s  In  these  local 
transfer  rates  and  consider  an  average  value  for  the  entire  sphere. 

The  tv  til  drop  In  either  temperature  or  partial  pressure  of  diffusing 
vapor  Is  assumed  tc  take  place  across  a  stagnant  film.  Consequently, 
heat  or  maos  transfer  equations  ouch  qj  given  in  the  firot  two  oll.dco 
are  used  hut  for  a  finite  film  thickness. 

The  heat  and  mass  transfer  equations  for  forced  convection  eva¬ 
poration  are  shown  on  the  next  slide.  It  Is  seen  that  these  equations 
are  similar  to  those  for  evaporation  in  a  quiescent  atmosphere.  The 
only  difference  is  the  use  of  a  Nusselt  number  to  account  for  the  effect 
of  forced  convection  on  evaporation  rate.  For  negligible  forced  con¬ 
vection  the  Nusselt  numbers  reduce  to  a  limiting  lever  value  of  2  and 
these  equations  become  identical  with  those  previously  shown  for  eva¬ 
poration  In  a  quiescent  atmosphere.  It  is  noted  that  the  heat  transfer 
IS  Cf  th?  Sf^lSpl.?  £02*2  thit  •  rtf  r*# 

the  diffusing  vapor  on  the  calculated  evaporation  rate. 

First,  let  us  consider  the  question  of  heat  trarsfer  for  spheres. 
Theoretical  approaches  to  heat  transfer  for  spheres  ’iave  been  made  by 
Johnstone,  Plgford,  and  Chapin  (ref.  22);  by  Kudryashev  (ref.  23);  by 
Drake,  Sauer,  and  Schaaf  (ref.  24);  and  by  Tang,  Duncan,  and  Schveyn* 

(ref.  25).  Predicted  values  of  average  Nusselt  number  over  the  sphere 
for  beat  transfer  in  air  are  shown  In  the  next  figure. 

Drake  followed  tho  assumptions  of  Johnstone,  but  used  a  different 
method  for  obtaining  the  solution  to  the  Fcurler-Polsscn  equation.  For 
Reynolds  numbers  below  1000  the  predicted  values  are  above  those  of 
Johnstone.  In  this  analyels  radial  velocity  was  assumed  to  be  aero  and 
tangential  velocity  was  assumed  everywhere  to  equal  mean  stream  velocity. 
For  Reynolds  numbers  approaching  aero  this  assumed  model  gives  a  limit¬ 
ing  Nusselt  number  of  2  as  is  required  for  heat  transfer  in  quiescent 
surroundings.  Kudryashev  obtained  an  analytical  solution  based  on  the 
thermal  boundary  layer. 

Another  approach  to  the  problem  of  heat  transfer  rates  from  spheres 
is  the  Reynolds  analogy.  This  con  be  stilted  that,  wh^n  heat  and  momentum 
are  transferred  in  cerresponding  ways,  tho  Nusselt  number  ever  the  product 
of  Reynolds  and  Prandtl  numbers  equals  half  the  friction  factor  &ti  deter¬ 
mined  from  skin  friction. 

_  Fcr  systea  such  as  flat  plates  where  skin  friction  Is  readily  deter-- 
mined,  the  Reynolds  analogy  holds  fairly  well  for  Prandtl  .numbers  cl  tee  to 
unity,  ChJ'iton  and  Colburn  (ref,  26)  suggested  simple  mcdlf Jcatlons  to 
the  Reynolds  analogy  for  Sctoldt  and  Frandti  numbers  other  than  unity. 
These  are  given  by  the  familiar  i -factors. 

For  the  case  of  bluff  bodies  where  it  Is  difficult  to  separate  akin 
friction  from  total  drag,  there  has  b»«n  little  work  <1  ..nr  to  dater-aine  the 
applicability  of -  the  analogy.  However,  worn  ly  Shrrvrcd  'rof .  27)  on 
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•  cyl  irJ'iii  e*.d,  r.:f-  ro: :r.tly.  work  V;.-  Tang  .v.\i  r-v. rk-  :v  (r<f.  £•“>)  m 
spheres  indicate  that  th*  analogy  ii,  lie  f*:ri j  ».‘U  icv  bluff  t:d1:s  for 
Prandtl  numbers  '•lea*  to  unity.  Ik-  Ku?»*li  r;.-.b.nr  br-s-.-.l  ::r.  the  n-.ynuda 
analog}'  and  calculated  ty  Tang  Uibiso  pr:*i:ntcd  l:.  thia  slide.  It  !.i 
noted  that  th-  friction  fartir  in  tVfe. calc-.laticn  la  baaed  on  the  sum¬ 
mation  of  local  friction  ovtr  thr  istlvs  surface-  aid.  not  cm  the  re&uituut 
of  friction  fore -a  in  the  direction  :f  tfca  freo  stream. 

Mow  let  ua  consider  empirical  correlations  for  heat  transfer  for 
spheres.  Ar.  extensive  sur’.ey  cf  published  data  cn  heat  and  naao  transfer 
to  spheres  was  conducted  by  Williams  (ref.  £8)  ir.  1042.  The  suggested 
relation  for  Hudee.lt  number  for  heat  transfer" la  given  by  equation  (1)  of 
the  next  slide.  Jn  the  range  of  lower  Reynolds  mr.be.ra  (<1C00)  the  data 
showed  considerable  scatter.  This  range  of  Reynolds  numbers  was  investi¬ 
gated  recently  by  Tang  (ref.  25),  who  determined  heat  transfer  rates  for 
steel  spheres  in  air  streams .  The  conditions  covered  included  ophere 
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and  air  streams  from  8°  to  100^  F  and  Reynolds  numbers  from  50  to  2G00. 

The  correlation  obtained  in  this  investigation  ia  given  by  equation  (2). 

Ran*  and  Marshall  determined  evaporation  rates  of  suspended  water 
and  benzene  drops.  Reynolds  numbers  ranged  from  zero  to  approximately 
200,  drop  diameters:  from  600  to  1100  microns,  and  air  temperatures  from 
room  temperature  to  approximately  400°  F.  For  evaporating  water  drops, 
the  data  were  correlated  by  equation  (3) . 

Ingebo  (ref.  29)  determined  evaporation  rato3  of  nine  pure  liquids 
in  air  streams  having  temperatures  ranging  from  approximately  80°  to 
1000°  F.  A  liquid  drop  was  simulated  by  means  of  a  wetted  cork  sphere 
having  an  average  diameter  of  0.69  centimeter.  Reynolds  numbers  ranged 
from  1200  to  5700,  In  a  later  report  (ref.  30)  data  were  obtained  In  air 
for  static  pressures  from  approximately  l/2  to  2  atmospheres.  The  majority 
of  these  later  experiments  vere  conducted  at  room  temperature.  In  order 
to  extend  the  correlation,  nr.  additional  li.nl ted  set  of  data  were  ob¬ 
tained  for  streams  of  helium,  argon  and  carbon  dioxide.  The  final  correla¬ 
tion  obtained  by  Ingebo  (ref.  30)  is  giver,  by  equation  (4).  Here  g  is 
the  gravitational  constant,  1  Is  the-  Sean  free  molecular  path,  anl  c 
the  root  moan  square  molecular  velocity. 

The  next  slide  presents  a  cenpariia-n  of  the  vwU.ua  predict*!  and  ex¬ 
perimental  curves  for  heat  transfer  for  spheres  in  tlv.  Here  Kusoelt 
uuiultu-  is  plultfcd  against  Reynolds,  ttebere  ii-r.%  *0  to  LOoO.  in  order  to 
compare  the  Rons  and  Marchall  eomlatu.r.  with  the  other  curves  on  the 
figure,  the  Prondtl  latmbor  in  thin  correlation  U<  taken  at*  that  cf  Ary  nir. 
The  correlation  obtained  t.y  ingebo  .1.s  based,  n.”.< -ng  <v Lh r.-r  factors,  on  a 
Schmidt  number  which  varl-s  vlith  the  liquid  being  vaporized.  The  Ingebo 
correlation  was  used  in  this  pgurc  by  taking  tn«  Schmidt  number  iw.l  ther¬ 
mal  conductivity  ratio  cquu'i  v..  unity.  The  correlation  factor  £l/c2  was 
calculated  cn  the  bonis  of  uir.  The  predict’d  .curves  of- Drake  and  Tang  are 
given  In  color. 
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It  la  seen  that  the  data  of  Tang  ore  more  closely  predicted  by  the 
Reynolds  analogy  than  by  the  curve  of  Drake  for  the  higher  Reynolds  num¬ 
ber  range.  For  the  lover  Reynolds  number  range  there  is  good  agreement 
between  the  predicted  curve  of  Drake  and  the  experimental  curves  of  both 
Ranz  and  Ingebo.  The  predicted  curve  of  Kudryashev  is  not  shown  on  the 
figure.  However,  it  lies  below  and  approximately  parallel  to  the  ex¬ 
perimental  curve  of  Tang. 

\ 

How  let  us  consider  the  treatment  of  forced  convection  evaporation 
in  terms  of  the  mass  transfer  equation.  Little  theoretical  work  has 
been  done  for  forced  convection  evaporation  specifically  in  terms  of  the 
mass  transfer  equation.  However,  in  view  of  the  analogy  between  heat  and 
mass  transfer  for  lew  temperature  evaporation,  the  theoretical  treatments 
fer  heat  transfer  could  be  applied  to  the  mass  transfer  case. 

A  theoretical  analysis  ror  mass  transfer  lu  nphcico  «a»  given  by 
Johnstone  and  Klelnachmidt  (ref.  31),  They  assumed,  the  sphere  contacts 
a  cylindrical  tube  having  an  inner  diameter  equal  to  that  of  the  sphere 
and  a  thickness  equal  to  the  distance  a  diffusing  molecule  could  travel 
in  the  time  the  sphere  moves  one  diameter.  Increase  in  relative  velocity 
between  sphere  and  aurrou:.  ling  gases  decreases  available  diffusion  time  and 
consequently  the  thickness  of  the  Imaginary  tube.  All  molecules  contained 
within  the  tube  are  considered  to  be  absorbed. 

Use  of  the  diffusion  equation  in  correlations  of  drop  evaporation  data 
Is  limited  to  evaporation  in  relatively  low  temperature  surroundings  Bince 
this  equation  becomes  highly  sensitive  to  experimental  errors  In  drop  sur¬ 
face  temperature  at  the  higher  surface  vapor  pressures ,  Correlations  of 
Nusaelt  number  for  mass  transfer  are  shown  in  the  next  slide, 

w'lilltuua  i  oi  related  mass  transfer  data  by  Frees sling :  .Tohnntntr*  nryi 
'.."'llama  (ref.  32);  Powell  (ref.  33);  and  McTnness  (ref.  31).  The  re¬ 
commended  equations  In  terms  of  the  Nusselt  number  for  mass  trunsfer  are 
given  by  equation  (l), 

Mai  a  el  and  Sherwood  (ref.  31 )  reported,  data  on  evaporation  of  water  und  • 
benzene  from  porous  spheres.  The  data  for  water  wire  In  good  agreement  with 
equation  (l).  Deviation  of  data  for  benzene  from  Mils  equation  was  thought 
to  be  attributed  to  incomplete  wetting  of  the  oj>h  re, 

Froeesling  correlated  data  for  room  temperature  evaporation  of  water, 
aniline,  naphthalene  and  nitrobenzene  for  Reynold h  numbers  from  4  to  ToZ . 
The  Froessllng  correlation  lo  given  by  ritual) en  (p). 

The  Ranz  and  Marshall  correlation  for  maun  transfer  in  given  ly  equation 
(3)..  This  correlation  is  of  same  form  an  their  heat  tmiinfer  ccrrelat.'.on, 
but  where  Prandtl  number  Is  replaced  by  llnlunidt  mini  •  i  .m  would  be  expected 
from  the  analogy  between  heat  and  maos  troiuitvr.  The  data  of  Frcedsling; 
Powell;  and  Maisel  and  Sherwood  agree  well  with  this  correlation. 
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Kinzer  and  C-unn  (ref.  2C)  obtained  a  predict-:-!  Nueselt  number  for 
Base  transfer  given  by  equation  (4).  The  equation  was  applied  to  eva¬ 
poration  rate  data  obtained  for  water  drop£  falling  freely  at  their 
terminal  velocities  under  rooa  air  conditions  Data  were  obtained  for 
drop  diameters  from  apcrcxlmately  20  to  4200  microns.  Application  of 
the  experimental  data  to  equation  (4)  Indicated  that  the  tern  F  In 
the  equation  varied  from  zero  to  a  value  greater  than  2  in  the  lev 
range  of  Reynolds  numbers .  For  Reynolds  numbers  from  approximately 
100  to  2000,  F  had  values  close  to  unity. 

It  la  seen  that  drop  evaporation  rates  for  forced  convection  have 
been  correlated  in  terms  of  either  the  heat  transfer  or  mass  transfer 
equation.  Use  of  either  of  these  equations  requiree  accurate  knowledge 
cf  drop  surface  temperature.  This  temperature  la  given  by  the  vet  bulb 
temperature.  For  the  special  ca3e  of  water  evaporation,  drop  surface 
temperature  may  be  taken  as  the  adiabatic  saturation  temperature  with 
little  error.  However,  for  other  liquids,  this  assumption  would  result 
in  serious  error  and  wet  bulb  temperature  must  be  used.  Vet  bulb  tem¬ 
peratures  may  be  calculated  by  simultaneous  solution  of  the  heat  and 
mass  transfer  equations.  For  the  case  of  evaporation  of  pure  liquids 
in  h::gh  temperature  surroundings,  Ingebo  has  obtained  a  correlation- for 
wet  bulb  temperature  in  terms  of  the  liquid  boiling  point. 

The  majority  of  the  data  have  been  obtained  for  pure  liquids. 

However,  Ranz  and  Marshal  1  determined  drop  temperatures  for  water  drops 
containing  dissolved  and  suspended  solids.  For  drops  containing  solids 
In  solution,  the  Initial  evaporation  rates  were  these  to  be  expected 
for  saturated  solutions  even  though  the  average  concentrations  In  the 
drops  were  belcw  the  saturation  value.  For  drops  containing  solids  In 
suspension,  the  Initial  drop  evaporation  rate  was  that  corresponding  to 
pure  water . 

There  appear  to  be  very  little  data  published  on  evaporation  cf  vide 
.boiling  range-  hydrocarbons.  However,  work  la  being  conducted  by  Iamb  and 
co-workers  In  this  field. 

The  final  case  of  single  drop  evaporation  to  be  considered  is  that 
of  unsteady  Btate  evaporation.  Here  the  Initial  drop  temper. t.irc  is  either 
above  or  below  Its  equilibrium  value  and  -  nr-equontly  ehangf-o  with  time. 
Both  theoretical,  and  experimental  work  has  teen  reported  for  this  case. 

The  analytical  treatments  generally  assume  either  a  finite  thermal 
conductivity  within  the  drop  as  given  by  that  of  the  liquid  or  asou-ue  an 
Infinite  thermal  conductivity  Circulation  vlthln  the  drop  rouults  J.n  an 
effective  thermal  conductivity  between  these  two  vhlurs,  the  circulation 
effects  being  more  pronounced  for  the  larger  drop  sizeu  and  Reynolds  rum-- 
bere  (ref.  37). 

Toppa  (ref .  38)  studied  evaporation  rates  cf  small  drops  falling 
through  a  high  temperature  atmosphere.  Initial  drop  diameters  ranged  from 
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approximately  300  to  550  microns ,  Hia  results  for  this  rather  narrow 
range  of  drop  diameters  indicated  that  the  evaporation  rate  varied  ap¬ 
proximately  as  the  drop  diameter  to  the  four  and  a  half  power  This 
is  contrasted  to  the  variation  of  evaporation  rate  vith  drop  diameter 
to  the  two  and  a  half  power  for  steady  state  evaporation  at  the  high 
Reynolds  numbers .  This  large  difference  was  attributed  by  Topps  to 
unsteady  state  conduction  of  heat  to  the  drop  interior  which  reduced 
the  heat  available  for  evaporation  of  the  smaller  drops, 

Godsave  (ref.  39)  investigated  theoretically  the  unsteady  state 
heating  of  burning  drops  of  heavy  fuel  oil-  -He  was  primarily  interes¬ 
ted  in  effects  of  the  unsteady  State  period  on  the  formation  of  carbon 
residue.  A  finite  thermal  conductivity  .was  assumed.  For  evaporation 
in  high  temperature  atmospheres  he  concluded  that  the  drop  interior 
tsii'.cl  to  remain  below  cracking '.temperature  .till  the  last  stages  of 
evaporation.  For  evaporation  in  lower  temperature  surroundingo  the 
drop  interior  tends  ito  reach  high  temperatures  in  the  early  stages  of 
evaporation  with  resultant  increases  in  cracking  and  cenosphere  forma¬ 
tion. 

Kinzer  and  Gunn  (ref.  3G)  used  an  ingenious  apparatus  to  determine 
the  time  rate  of  change  in  the  average  temperature  of  falling  water  drops. 

A  simplified  theoretical  expression  was  obtained  for  the  time  required 
for  the  drops  to  reach  63  percent  of  the  equilibrium  temperature  difference. 
The  analysis  assumed  thermal  conductivity  of  the  drop  interior  to  be  in¬ 
finite  .  Satisfactory  agreement  was  found  between  the  predicted  and  ex¬ 
perimental  value. 

Recently  El  Vakil,  Myers,,  and  Uyehara  (ref,  40)  studied  the  unsteady 
state  evaporation  of  pure  liquid  drops  for  conditions  similar  to  those 
existing  in  Jet  engine  combustors.  'In  their  analysis,  thermal  conductivity 
of  the  drop  interior  was  also  assumed  to  be  infinite  Temperature -time, 
mass-time,  and  penetration-time  histories  were  calculated  for  drops  eva¬ 
porating  in  air.  This  type  of  calculation  involves  a  step  by  step  trial 
and  error  procedure  for  simultaneous  solution  of  the  heat  transfer,  maos 
transfer,  and  drop  motion  equations.. 

A  typical  figure  is  shown  in  the  next  slide.  Here  percent  evaporated, 
relative  velocity,  and  drop  temperature  are  plotted  against  time  in  milli¬ 
seconds.  The  calculations  are  for  an  iBoor.tnne  drop  having  an  initial  dia¬ 
meter  of  50  microns,  initial  temperature  of  no0  p.  and  initial  vei ori'ty 
relative  to  air  of  100  feet  per  second.  The  air  was  assumed  to  be  at  a 
pressure  of  l/2  atmosphere  and  temperature  of  1000°  F.  For  these  condi¬ 
tions,  the  amount  of  heat  going  to  the  drop  Interior  far  exceeded  that 
supplying  latent  heat  of  vaporization  during  the  major  portion  of  the 
unsteady  atate  period.  Calculations'  of  this  type  Indicate  that  at  high 
air  temperatures  and  for  high  volatility  fuels,  the  unsteady  state  period 
nay  represent  an  appreciable  part  of  the  total  vaporization  time 

We  now  come  to  the  vaporization  of  liquid  sprays.  Since  the  subject 
of  oproy  drying  will  be  treated  in  another  paper,  thi3  portion  will  be 
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limited  to  a  brief  discussion  cf  data  of  mere  direct  concern  tc  the  com¬ 
bustion  chamber  designer. 

Extension  of  single  drop  data  to  sprays  ia  difficult  since  both  drop 
size  distribution  of  the  spray  and  relative  velocities  between  the  air  and 
droplets  must  be  known.  In  addition  there  are  such  complicating  factors 
as  drop  distortion,  unsteady  state  evaporation,  and  interaction  between 
drops. 

Several  theoretical  analyses  of  spray  evaporation  have  been  made  by 
ueing  single  drop  relations  in  conjunction  with  assumed  drop  aize  distri¬ 
butions.  One  such  analysis  is  that  described  by  Probert  (ref.  41).  In 
this  treatment,  all  drops  were  assumed  to  have  zero  relative  velocity  with 
respect  to  the  air.  The  spray  was  assumed  to  follow  the  Rosin-Rammler 
distribution  onn«+.inn.  with  this  analysis,  the  fraction  of  unevanorated 
spray  after  a  given  period  of  time  is  given  in  terms  of  the  opray  drop 
size  distribution  constants  and  the  evaporation  rate  for  a  Bingle  drop  in 
a  quiescent  atmosphere. 

Another  theoretical  treatment  of  spray  vaporization  has  been  presented 
by  Tribus,  Klein,  and  Rembowskl  (ref,  42).  Zero  relative  velooity  between 
the  air  and  droplets  was  also  assumed.  It  is  shown  that,  given  an  initial 
drop-size  spectrum,  successive  spectra  may  be  deduced  by  calculating  the 
variation  in  drop  size  of  the  largest  drop  and  relating  all  other  drop 
sizes  to  this  maximum  drop  size. 

Because  of  the  lack  of  data  on  drop  size  distribution  for  a  particular 
apparatus,  .it  is  generally  difficult  to  apply  evaporation  rate  relation¬ 
ships  for  single  drop  to  sprays.  Some  attempts  have  been  made  to  obtain 
direct  measurements  of  spray  vaporization  rates  in  air  streams.  An  in¬ 
vestigation  of  this  type  has  been  described  by  Fledderman  and  Hanson  (ref. 
43).  The  Influence  of  turbulence  and  air  velocity  on  the  spray  vaporization 
rate  was  studied.  A  theoretical  analysis  of  spray  vaporization  was  also 
Included  in  this  study,  In  the  analysis,  Froessling's  equation  for  single 
droplet  evaporation  was  combined  with  the  drop  size  distribution  equation 
of  Nukiyoma  and  Tanasawa  (ref.  44). 

The  experimental  measurements  were  obtained  by  sampling  hexane  sprays 
at  various  distances  downstream  of  a  hollcw-cone  spray  nozzle.  The  percent 
of  evaporated  spray  was  determined  over  a  limited  range  of  air  velocities. 
The  results  indicated  a  strong  dependence  of  evaporation  rate  on  the.  rela¬ 
tive  velpclty.  Increased  turbulence  intensity  was  also  found  to  Increase 
the  spray  vaporization  rate.  The  test  rcsulto,  however,  did  not. provide 
conclusive  Information  on  the  effect  of  scale  of  turbulence  on  tho  evapora¬ 
tion  rate. 

The  results  of  a  few  measurements  cf  the  evaporation  rate  of  kerosene 
Bp?3y  3  la  ecccat1.cLi.ly  ctlll  air  hare  bcea  reported  ty  Sticks  ( ref  *  45)* 


Ia  this  program,  the  kerocene  was  Injected  from  a  swirl  chamber  fuel 
nozzle  at  pressures  of  50  and  00  pal.  Over  this  range  of  oondltlona, 
from  0.2  to  0.3  percent  of  the  a  pray  vaa  found  to  evaporate  per  second. 

This  la  contrasted  to  a  calculated  value  or  approximately  50  percent 
per  second.  These  predicted  values  vere  obtained  by  using  the  Langnulr 
equation,  Probert's  relations,  and  drop  size  distribution  values  ob¬ 
tained  from  the  Boven  and  Joyce  relations  (ref.  4G).  This  very  large 
difference  betveen  predicted  and  experimental  values  vaa  attributed 
primarily  to  the  Inability  of  the  Langmuir  equation  to  predict  eva¬ 
poration  rates  for  single  drops  when  in  a  cloud  of  drops. 

A  study~of  the  effect  of  the  air  flow  and  fuel  injection  parameters 
on  the  evaporation  of  gasoline -type  fuel  sprays  has  recently  been  made 
by  Bahr  (ref.  47).  Thla  investigation  vaa  conducted  over  ranges  of  con¬ 
ditions  common  in  raa-Jet  and  afterburner  engines.  An  empirical  cor¬ 
relation  vaa  obtained  ror  tne  evaporation  over  the  range  of  conditions 
Investigated.  This  study  will  be  discussed  by  Mr.  Bahr  at  the  conclusion 
of  this  paper. 

In  summary,  ve  find  good  agreement  betveen  predicted  and  experimental 
steady  state  evaporation  rates  of  single  drops  in  low  temperature  quies¬ 
cent  surroundings.  Experimental  data  on  evaporation  in  high  temperature 
quiescent  surroundings  are  limited  to  the  case  of  burning  drops.  Again 
good  agreement  is  found  betveen  experiment  and  theory. 

Theoretical  approaches,  to  steady  state  evaporation  under  forced  convec¬ 
tion  are  far  more  difficult'.  Consequently,  a  number  of  simplifying  assump¬ 
tions  have  been  required  in  order  to  obtain  solutions.  Agreement  betveen  ' 
theory  and  experiment  is  satisfactory  in  viev  cf  the  approximations  used 
in  the  theory.  Semi -empirical  correlations  for  steady  state  forced  convec¬ 
tion  evaporation  arc  available  in  terms  cf  either  the  jjgivf  woo  n  trs!isf?r 

equations.  Drop  eurface  temperatures  may  be  determined  from  simultaneous 
solution  of  these  equations  or  obtained  from  a  correlation  in  terms  of  the 
liquid  belling  point  for  the  higher  temperature  evaporation  conditions. 

The  major  portion  of  the  experimental  data  vaa  obtained  for  pure 
liquids.  Some  data  vere  obtained  cn  the  effect  of  both  dissolved  and  sus¬ 
pended  solids  on  the  evaporation  rate  cf  vater  dropB.  There  appear  to  be 
very  little  data  available  at  the  present  time  that  vould  permit  predic¬ 
tions  to  be  made  for  .he  evaporation  rate  of  vide  boiling  range  hydro¬ 
carbons  . 

Calculations  have  been  made  for  the  unsteady  state  period  of  drop  eva¬ 
poration.  Seme  of  these  calculations  indicate  that  unsteady  state  evapora¬ 
tion  may  represent  an  appreciable  portion  of  the  total  drop  evaporation 
time.  However,  there  ore  very  little  experimental  data  available  for  this 
part  of  the  vaporization  process. 


Several  theoretical  approaches  have  been  made  to  the  problem  of 
predicting  evaporation  rates  of  liquid  sprays.  However,  hut  a  limited 
amount  of  experimental  work  has  been  dene  to  determine  how  well  data 
for  single  drops  would  apply  to  vaporization  rates  of  sprays  in  air 
stream.  This  type  of  investigation  would  involve  the  difficult  task 
of  measuring  simultaneously  both  drop  size  distribution  and  evaporation 
rate.  One  investigation  indicated  an  appreciable  difference  between 
predicted  and  experimental  spray  evaporation  rates  in  a  low  tempera¬ 
ture  low  velocity  air  stream.  An  empirical  correlation  was  obtained 
from  direct  measurements  of  evaporation  of  an  isooctane  spray  in  a  high 
velocity  air  stream.  Further  experimental  work  of  this  type  is  needed 
to  determine  the  effect  of  fuel  properties  and  fuel  injection  systems 
on  spray  evaporation  rates. 
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EVAPORATION  OF  SINGLE  DROPS' 

1.  QUIESCENT  SURROUNDINGS 

(a)  LOW  TEMPERATURE 

(b)  HIGH  TEMPERATURE  , 

2.  FORCED  CONVECTION 

( a )  HEAT  TRANSFER 
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3.  UNSTEADY  STATE 
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EFFECT  OF  ENTHALPY  CHANGE  OF  VAPOR  ON  PREDICTED 
EVAPORATION  RATES  OF  WATER  AND  ISOOCTANE  DROPS 
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COMPARISON  OF  PREDICTED  MJSSELT  NUMBERS 
HEAT  TRANSFER  FOR  SPHERES 
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EUPIRICAL  CORRELATION  OF  MEAi-T^/.NIir-T  DATA  FOR  SPHERES 
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EMPIRICAL  CORRELATIONS  FOR  MASS  TRANSFER  FOR  SPHERES 
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INJECTION  HI  AFTERBURNING  TURBOJETS 


I  Introduction 

From  the  very  beginning  of  Jet  Propulsion  the  operational  limitations 
of  the  simple  turbojet  were  understood  and  discussed,  and  methods  for 
improvement  were  suggested  (Ref.  l). 

Among  others,  these  limitations  in  performance  included  the  following: 

a*  Relatively  long  take-off  runs} 

b.  Low  initial  rates  of  climb} 

c.  The  requirement  for  bettor  top  speed  performance  (short  duration). 

Afterburning,  as  it  was  first  introduced  (Ref.  l),  also  included  cooled 
compression  in  the  compressor  (e.g.  by  liquid  injection  In  the  inlet).  This 
increases  the  thrust  of  the  unit  by  increasing  the  temperature  of  the  gases 
before  expansion  in  the  propelling  nozzle  and  by  decreasing  the  compression  work 
In  the  compressor.  Figure  1  shews  this  cycle.  Considerable  theoretical  and 
experimental  work  was  then  performed  in  order  to  realize  the  above  advantages* 

It  should  be  stressed  here  that,  although  considerable  improvements  C3n 
bo  obtained  in  the  subsonic  flight  region  by  afterburning  (Ref.  2),  its 
employment  at  supersonic  flight  speeds  has  been  found  to  be  not  only  desirable 
but  imperative  (Ref.  J). 

Although  the  idea  of  afterburning  is  simp..*,.  its  full  realization  in  an 
engine  is  fraught  with  a  number  of  difficulties  each  asv 

a.  Stable  operation  at  any  altitude  without  screech} 

b.  Operation  at  high  altitude  without  flame  out} 

e.  Propelling  nozzle  cooling. 

The  fir  at  two  difficulties  are  strongly  influenced  by  the  injection 
process  of  the  fuel  in  the  afterburner,  ard  will  be  discussed  in  more  detail. 

Generally  speaking,  two  extreme,  although  representative,  operational 
eases  of  afterburners  are  apparent: 

1.  The  afterburner  operates  only  a  comparatively  short  part 
of  the  total  flight  time} 

2.  The  afterburner  operates  almoot  the  whole  flight  time. 
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The  operational  requirements  of  the  first  case  are  as  follows i 

iv^Short  take-off  run; 

b.  High  pate  of  climb; 

c«  Short  tina  -of  top  thrust  operation) 

d,  long  time  economical  cruising. 


As  the  afterburning  system  will  operate  only  during  a  short  time 
Interval  it  should  be  designed  in  such  a  way  as  tsTavoid  large  losses  duri "g 
the  long  non-afterburning  time*  In  this  case  higher  afterburner  losses  cay  bo 
tolerated  during  the  afterburning  period* 

Tho  operational  requiremsnts  of  the  second  case  are  es  follows! 

a*  Long  time  economical  cruising  with  afterburning  ct  high  Hich  speeds 
and  altitudes; 

b*  Good  take-off  and  transonic  operation  is  desirable; 

e*  The  operational  altitude  and  speeds  are  well  above  those  of  the 

first  ease* 


Rrom  the  above  it  is  obvious  that  the  afterburning  aystea  should  be 
designed,  so  as  to  operate  with  minimum  loss  during  cruising.  In  this  case  higher 
afterburner  losses  may  be  tolerated  during  the  no  a- aft  or  burning  period. 

The  high  altitudes  and  Uach  speeds  anticipated  may  put  severe  restrictions 

on  the  design. 

Thfl  n<m®  stability  at  cruising,  and  during  the  acceleration  phases 
should  have  an  amply  safe  range  of  operation;  herrovor,  that  may  not  be  always 
feasible  and  fl.arae  out  at  altitude  may  bacons  a  serious  problem. 


II  Types  of  Afterburning 


It  la  well  known  (Ref.  h)  that  the  combustion  of  a  fuol  in  a  propulsion 
system  can  be  effected  according  to  the  following  modest  V  ^  ° 


«*  Injection  of  liquid  fuel  in  a  gas  otream; 
b*  Injection  of  vaporised  fuel  in  a  gas  stream; 
c.  Catalytic  combustion  on  a  suitable  hot  surface. 


We  shall  discuss  here  only  the  mode  of  liquid  injection*  Afterburning 
types  are  sub-divided  according  to  the  position  of  the  injection  apparatus  and 
the  direct  ion -of  injection  relative  to  the  gas  stream* 

According  to  the  above,  the  following  main  types  of  afterburning  may 
distinguish  themselves t 

a*  Preturbine  injection*  This  type  has  been  investigated 
experimentally  by  the  Air  Farce  (Power  Plant  Laboratory, 

Weight  Patterson  Air  Force  Base)  (Fig*  2-a.)j 

b*  After  turbine  upstream  injection  with  flaae  holder 
stabilization  (Fig,  2-b) j 

e*  After  turbine  upstream  injection  with  Haae  in  the 
trailing  vortices  of  the  turbine  blades  (Fig*  2-c»)j 

d*  After  turbine  downstream  injection  with  flame  holder 
stabilization  (Fig*  2-d)j 

a*  After  turblio  downstream  injection  with  tail-cone-end 
injector  (Fig.  2-e.)«, 

i 

XIX  Advantages  and  Disadvantages  of  the  Different  Typest 

The  first  type,  namely,  preturbine  injection,  does  not  need  a  flame 
holder  'as  the  flame  stabilization  is  mostly  effected  by  the  turbine  blades* 

This  type  when  properly  designed,  is  accompanied  by  high  oombustion  efficiency 
and  a  short  length  of  propelling  nozzle  is  always  attained.  The  main 
disadvantage  of  this  typo  is  complete  lack  of  control  of  the  temperature 
distribution,  which  results,  in  the  appearance  of  hot  spots  but  no  combustion 
instabilities.  Cooling  of  the  fuel  injectors  may  be  necessary. 

The  second  type  of  afterburner,  when  properly  designed,  gives  high 
combustion  efficiency  and  a  satisfactory  temperature  distribution.  A  short 
length  of  propelling  nozzle  is  always  obtained.  However,  because  of  the 
presence  of  flame  holders,  choking  of  the  flow  in  the  flame  holder  region 
appears  at  high  fUel-air  ratios,  accompanied  by  the  deleterious  oscillativig 
phenomenon  known  as  screech* 

The  third  type  seems  to  possess  the  best  all-around  characteristics, 
namely,  high  canbustlon  efficiency  and  controllable  temperature  distribution. 
Since  there  are  no  flame  holders,  no  flow  blockage  exists  and  consequently, 
it  Is  quite  free  of  screech*  This  typo  may  need  considerable  cooling  of  the 
fuel  injectors*  The  length  of  the  propelling  nozzle  can  be  kept  short* 

•The  fourth  type  requires  a  comparatively  longer  length  of  afterburner 
to  keep  the  combustion  efficiency  at  a  permissible  level.  This  sometimes  is 
accompanied  by'high  pressure  leases.  Because  flame  holders  are  employed  for 
flame  stabilization,  blocking  of  the  flow  occurs  and  this  type  is  very 
susceptible  to  screech* 
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The  fifth  type  is  the  oldest  utilized  experimentally.  The 
combustion  efficiency  obtained  in  19^5  around  60 %,  Alter  the  war,  by 
thb  use  of  long  propelling  nozzles,  the  combustion  efficiency  was  improved 
at  the  expense  of  high  pressure  losses*  This  type,  when  properly  designed, 
is  not  susceptible  to  screech* 


It  is  evident  that  the  shorter  the  length  of  the  propelling  nozzle 
the  easier  the  solution  of  the  cooling  problea* 


Table  I  gives  a  summary  of  the  above  discussions 
TABLE  1 


Injection 

Type 

Temperature 

Distribution 

Combustion 

Efficiency 

Pressure 

Losses 

Screech* 

Propelling 

Nozzle 

•  1 

Bad 

Good 

Low 

Free 

Short 

2 

Good 

:  Oood 

Lew 

Badly 

Short 

3 

Oood 

.'Oood  , 

Low 

Free 

Short 

U 

toy  Be  Oood 

» 

Hot  So  Oood 

High 

Very 

Badly 

Long 

5 

toy  Be  Good 

Poor 

Hot  So 

High 

May  Be 
Free 

Very  Long 

In  the  case  where  a  part  load  operation  is  required,  staggering  of 
the  injectors  Into  two  or  three  rnri  results  in  a  satisfactory  operation  of  the 
afterburner* 

IT  Injection,  Evaporation  and  Ignition 

The  main  objective  of  Injection  and  ignition  is  to  occurs  a  unifbrn 
temperature  and  velocity  distribution  during  and  after  eoobustion  and  also  to 
give  sufficient  roaldonce  time  to  the  fool  droplets  so  that  the  process  of 
evaporation,  airing  and  Ignition  will  have  enough  time  to  be  effectively 
performed. 

Sr  order  to  provido  a  better  design,  quantitative  data  on  the 
injection,  evaporation,  nixing,  and  ignition  of  fuel  droplets  are  necessary. 
Ignition  data  are  mostly  concerned  with  spontaneous  ignition. 

The  tine  end  length  of  droplet  penetration  are  Important  and  can  be 
calculated  iron  taoom  formulas  both  for  single  droplets  and  for  simple, 
mathematical  droplet  distributions  (sprayrj)* 


Upstream  injection  is  not  as  sensitive  to  droplet  size  and 
distribution  as  is  downstream  injection.  This  results  from  the  larger 
droplet  having  greater  upstream  penetration.  Again  greater  penetration  is 
accompanied  by  longer  residence  tine.  The  opposite  is  true  for  det/nstream 
injection.  Under  similar  flow  conditions  such  as  velocity,  flame  holder 
distribution,  temperature,  etc,,  the  deterioration  of  altitude  operation  of 
an  afterburner  with  downstream  injection  is  more  pronounced  than  that  of  one 
with  upstream  injection  because  of  shorter  residence  time* 

The  variables  appearing  in  the  design  of  an  afterburner  are: 

a.  Temperature  of  the  turbine  outlet  gases) 

.  b.  Pressure  of  the  turbine  outlet  gases) 

c.  Oxygen  depletion  (vitiation)  of  the  turbine  outlet  gases) 

d.  Fuel  characteristics) 

; e.  Fuel  pressure) 

f.  Fuel  droplet  size  and  distribution. 

Work  performed  by  B.P,  Hull ins  during  uhe  last  decade  at  the  N.G.T.E. 
indicates  that  the  influence  of  the  droplet  size  and  the  fliel  characteristics 
(when  kerosene  type  fuels  are  Used),  do  not  have  a  pronounced  influence  on  the 
Ignition  lag.  All  other  variables  have  a  profound  effect  on  the  evaporation  and 
ignition  process. 

Figure  3  shows  the  ignition  lag  of  a  standard  fuel  plotted  as  a 
function  of  gas  pressure  at  constant  gas  temperatures.  It  is  cbvious  that  as 
the  pressure  decreases  for  a  given  gw  temperature,  the  ignition  lsg  Increases 
rapidly,  thus  resulting  in  a  rapid  decrease  of  combustion  efficiency  with 
altitude.  This  may  be  remedied  up  to  a  certain  degree  by  an  increase  of  the 
turbine  outlet  gas  temperature  or  an  increase  of  the  flame  holder  drag. 

The  above  indicates  quite  clearly  a  ceiling  in  afterburner  operation 
upon  which  the  injection  type  has  a  profound  effect.  In  spite  of  the  above,  it 
is  encouraging  that  the  introduction  of  turbine  blade  cooling  results,  to  some 
extent,  in  higher  turbine  outlet  temperatures  and  thus  raises  the  operating 
ceiling  of  the  afterburner  by  reducing  the  ignition  lag. 

Figure  U  shows  the  influence  of  vitiation  (contamination  with 
combustion  gases).  2h  this  the  ignition  lag  is  plotted  against  the  relative 
depletion  of  oxygen  at  constant  gas  temperatures.  This  figure  indicates  that 
when  maintaining  all  other  conditions  constant,  increasing  the  turbine  outlet 
teoyerature  will  result  in  a  higher  degree  of  vitiation  which  will  tend  to 
decrease,  somewhat,  the  advantages,  shown  in  Figure  3,  of  an  increase  in 
temperature. 
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7  Beat  Addition;  Stability;  Screech 

Generally  speaking,  afterburning  may  be  considered  as  a  heat 
addition  process  in  a  more  or  less  high-speed  gas  stream.  The  one  dimensional, 
steady-flow  theory  of  heat  addition  in  a  variable  area  duct  with  losses  has 
been  given  before  (Ref.  5  and  Ref.  6).  The  conclusions  of  this  theory  may  ba 
expressed  as  follows  t  Heat  can  bo  added  in  a  duct  in  a  stead/  fashion*  or  in 
the  more  realistic' case  in  a  stream  tube,  until  the  local  Bach  number  of  the 
flow  teaches  unity.  Then  the  total  head  pressure  loss  and  the  total  temperature 
reach  an  extreme  value.  As  it  was  shown  in  the  above-mentioned  reference,  the 
variation  of  total  head  pressure  loss  is  as  follavsi 


+  #,y 

D  o  I  Ti  Si.00  A  J  '  ' 


Where 

fb  -  total  pressure 

Tt=  total  temperature 

?  s  specific  heat  ratio 

M  -  llach  number 

A  = cross-sectional  area 

C^,=  friction  coefficient 

@s slope  of  the  meridian  of  the  duct 


This  formula  suggests  that  in  order  to  add  heat  under  low  pressure  losses,  the 
local  Mach  number  of  the  flow  should  be  kept  lew.  This  can  only  be  realized  by 
an  increase  of  the  area  of  the  3tream  tube  or  the  duct  area.  Such  an  increase 
may  have  an  adverse  effect  on  the'  frontal  area  of  the  propulsion  system  and 
afterburner  cooling. 

The  main  assumptions  of  the  above  developed  theory  were  that  the  How 
is  steady  and  the  velocity,  temperature,  etc.,  distributions  were  nearly  uniform, 
thus  enabling  us  to  treat  the  problem  as  one  dimensional.  Experience  shows, 
however,  that  neither  of  those  two  assumptions  i3  realistic  near  the  sonic  velocity 
because  present  injector  designs  do  not  yield  uniform  distributions  across  the 
duct. 
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It  is  well  known  that  at  the  relatively  high  temperatures 
encountered  in  afterburners,  the  combustion  is  unsteady  and  oscillatory,  the 
combustion  taking  place  air.  ys  in  the  turbulent  region. 

Again  because  of  non-uniform  distributions  of  velocity  and 
temperature,  especially  near  solid  obstructions  such  as  flame  holders,  the  Mach 
number  distribution  is  non-uniform  and  consequently  at  certain  stream  tubes 
where  fuel  is  available  for  combustion  the  Mach  number  may  approach  unity. 

It  is  well  known  (Ref.  7)  that  no  steady  transition  can  occur  from  a 
deflagration  (subsonic  energy  release)  to  a  detonation  (suporeonic  energy  release) 
the  transition  always  occurs  in  an  oscillatory  fasHon.  This  transonic  energy- 
release  has  been  called-spinning  detonation  and  has  been  treated  theoretically 
by  Uanson  (Rof.  6).  It  has  been  found  that  spinning  detonation  is  closely 
connected  with  "screech" • 


Dttscrapuxuu  tuu  piiyuxuax  t*Apx«ukivxuu  ux  uto  i/runauiixC  OaOxxxaw*/ 

phenomenon  known  as  "screech"  was  given  by  this  author  in  January  1950.  Similar 
explanations  with  more  mathematical  details  were  given  later  by  Vfitbeck  and 
by  Dennis  and  Putnam. 

"Screech"  is  a  transversal,  transonic  oscillation  which  results 
from  non-uniform  transversal  distribution  of  velocity  and  Mach  number.  This  non- 
uniformity  of  velocity  and  temperature  is  produced  and  accentuated  by  the 
injection  and  combustion  processes.  Consequently,  the  deleterious  phenomenon  of 
screech  oan  be,  avoided  or  delayed  by  efficient  design  of  the  afterburner  In 
such  a  way  that  the  injection  and  combustion  processes  will  produce  a  uniform 
velocity  and  temperature  field. 

VI  Matching  and  Acceleration 

^hca  variable  degrees  of  afterburning  are  used  and  the  temperature 
of  afterburning  and  the  area  of  propelling  nozzle  are  known,  the  "equilibrium 
lines"  of  operation  may  be  found  if  only  one  such  line  is  known. 

As  is  shown  in  the  Appendix,  the  "equilibrium  running"  point  on  the 
matching  diagram  (Figure  5)  corresponds  to  a  fictitious  non-afterburning  turbojet 
with  a  fictitious  temperature  ratio. 


(2) 


Firom  the  above,  wo  may  conclude  that  when  the  propelling  nozzle  area  k,  stays 
constant,  aftorburning(n3*/,T^  results  in  an  increase  of  the  fictitious 
temperature  ratio ('l\t/Tt^u3  raising  the  equilibrium  line  of  operation,  and,  if 
no  sufficient  margin  is  available,  surging  of  the  compressor  may  result.  This 
predicament  can  be  remedied  however,  by  an  increase  of  the  propelling  nozzle 
area,  as  is  shown  in  the  above  formula.  From  the  above,  it  is  obvious  that  after¬ 
burning  has  the  same  effect  as  closing  the  propelling  nozzle.  During  acceleration 
the  operation  point  moves  toward  higher  temperature  ratio,  values;  consequently, 

Seat  care  should  be  given  during' the  design  of  acceleration  controls  of  tho 
Jection  system  tc  keep  the  point  of  operation  within  tho  stability  limits. 


In  the  case  of  variable  area  propelling  nozzles  this  can  be  effected  by 
linking  the  acceleration  controls  with  the  propelling  area  controls. 

Til  Conclusions 

The  infection  process  has  a  great  influence  upon  the  operation 
of  afterburning  turbojets.  With  the  present  lack  of  fundamental  knowledge  on 
injection,  empirical  design  of  afterburners  is  followed  by  a  number  of  serious 
limitations,  such  as  altitude  flame  out,  screech,  etc. 


APPENDIX 


Equilibrium  running  characteristics  of  afterburning  turbojets. 

The  method  followed  was  indicated  in  Ref.  The  mass  gas  flow  equation 
through  the  propelling  nozzle  may  be  given  as  follows: 


Where: 


Qs“ 


a*Tk 

A*  Rt 


mass  flow  rate  through  the  afterburner 


TV  afterburner  maximum  total  gas  temperature 
afterburner  total  pressure 


propelling  nozzle  area 
prepolling  nozzle  cutlet 
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flrom  the  above,  the  follovrLng  may  be  obtained: 


Where:  CJ,  *  air  mass  flew  through  the  compressor 

Tit"  total  temperaturo  at  the  compressor  inlet 

pressure  at  the  compressor  inlet  and  outlet 
total  pressure  before  the  turbine 
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•Then  the  total  pressure  loss  in  the  main  combustor  and  afterburner  are 
the  following  equations  hold  true: 


Fk  —  Hi  ±  .Si  = 

Rt  -Rt  P* 

From  these,  solving  for  Bt  ire  obtain: 


JRl  _  _  ASi  —  .  Ot .  I  Tst .  .Sa  .  Si  _ _ ! - - 

P,t "  Rt  Rt  "  P?t  Q,  P«  R* 

In  the  case  of  turbine  choking  or  near  choking  the  following  equation  holla 
true: 

CU  . .  /all' 

Ptt"“'/'n, 


(3i 


(U) 


Where: 


T*iT*-  the  inlet  and  outlet  turbine  total 
1  temperature  respectively 

F^t  V?t«  the  iiilet  ani  outlet  turbine  total 
pressure  respectively 


Consequently,  we  find: 

PU-Q./tS  Qs  pL&.[y>2l\  4  ^R«\ _ K _ !_+  ABa 

Rt"  Rt  V  Rt)  Rt. 

In  a  single  non-afterburning  turbojet  the  fc llc.ring  holds  truo: 


P,t  .Q/iP  C2a„® 

Rt  “  Rt  Qi  i  T,i.  KfK) 


r\t 


(7) 


Cy  cosparing  the  above  equations  it  may  be  seen  that  the  same  characteristic 

m-,..  =aj  ....  voJw  r/.w  «■ 

ov^ucouuh  uvxu^  iv*  vtit<  ui  vv 

is  considered,  as  follows: 


-  - .  <V« m  t  Va  a /v  a *»Vt t «>* m 4  m  ••  t « i v»V\ a 4  4  r  it  f  4  a+  4  f  *  Aito  4  AniAnt* nf  ii  yss  >tof  4  n 

o^uaouuii  UVXU«;  iv*  vut*  uivut^ut  vtt*  wwjw  u  a4.w^v*vmw  vwMr>-.  UtUrg 


//3ST=  .  pS  ./ET.f, + AS»\f  -MMO 

y  Tit/  o,  rnt  fnt  v  pjt 

Whore  tho  superscript  aero  Indicates  non-afterburning  and 
^3  “  gas  mass  Hot  rate  through  the  turbine 


(8) 
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